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a^opsis 


For developing countries like India, where material 
to labour cost is high , shell foundation for columns will 
play a vital roLe in saving material m cons tiuc li on „ Shells 
on elastic foundations me of recent; origin end no realistic 
analytical work has been done as yet although, some experi- 
mental work has boon reported m tlio literature. A rovaow 
of the work done on sheila on elastic fo midations has been 
given m Chapter I, 

Various foundation models have boon proposed by 
several authors. A comprehensive review of the various 
foundation models has also been given m Chapter I, The 
earliest and simplest formulation of elastic foundation was 
by Winkler who assumed that the foundation model consists 
of closely spaced independent linear springs. This model 
lacks continuity. Moot cons .Is lout formulation would be to 
consider the foundation as a semi-inf in i to elastic continuum, 
but the analysis in this case is very comp Lex. Several 
authors have pro posed foundation models which arc a com- 
promise between the above mentioned two extremes. These 
foundation models involve more than ono parameter to charac- 
terise the foundation medium. To mention a few we have 
models proposed by Wieghar d b , Filonenko-dorodich , Pasternak, 
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Reissner , Ilotenyi, Vlasov and leontov a o cl kameswar Rao , 

. Das and Anandaknshnan . 

Using variational method a foundation model m 
curvilinear orthogonal coordinates has been proposed m 
Chapter II, Thu foundation model includes displacements 
m all the three orthogonal coordinate directions, hence 
making it possible' to use the proposed model for foundation 
of any type of stricture and m any set o f orthogonal 
coordinate system. This model is capable of analysing 
foundations of both finite and infinite thicknesses. Multi- 
segment technique, which has been used to solve the boundary 
value problems has also boon described in this chapter. 

In Chapter III pro posed foundation model has been 
used to analyse thin shallow spherical and conical .shells 
on elastic foundations under axisymmo t rie normal load, 
transferred by a column. Three important types of boundary 
conditions namely simply supported, L'i.?;od or built m and 
free are considered. In case of free boundary condition 
the flank portion (foundation beyond thr shell ) has been 
taken into account. Depth of foundation has been considered 
as finite m tho numerical examples. In case of spherical 
shells, two parameters accounting for thinness and shallow- 
ness of the shell, are considered. In case of conical shell 
only one parameter which accounts for the shallowness of the 
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3 lib 11 has been cone. id© rod . .Results are ptesentod m non- 
dimensional form bo that they enn be urn d for design pum 
poses. Curves for moments, displacement t, mplone forces 
and foundation reactions art presented Lor different values 
of parameters. 

Ln Chapter IV, proposed loun.dalj.on model has been 
used to analyse tli lu shallou spherical and conical shells 
on clastic fouuda l.u ons subjected to onti symmolr ic loads, 
loads m radial dime Lion and moments io radial direction 
transferred by tlic column have boon considered. Simply 
supported, fixed and free boundary conditions have boon 
used. Depth of foundation has been tsbeu as finite in the 
numerical examples that arc solved. Some parameters as 
mentioned in Chapter III are accounted Joe. The results 
arc presented m nondunons lonai form. Ouives for moments, 
inplono forces, din placements and foundation rone tions 
have boon presented for different values of. par an e tors » 

Conclusions and suggestions for in lure research have 


been made in Chapter V. 



CHAPTER 1 


INTRODUCTION 


1.1 GENERAL: 

Por developing countries like India, where material 
to labour cost is high, shell foundation for columns will 
plays a vital rolo in saving material in construction works. 
Shell structure Cor roofs have become very popular for last 
25 years because of the inherent property of its being 
structurally more efficient than any other type of roof 
structure . 

Shells on elastic foundations is of recent origin. 
Shell foundation is of considerable importance in structural 
engineering because of many reasons. To mention a few 

I 

(.i) the material saving is upto 50 percent in many a casoa 
of construction as pointed out by Ruhlo [1] and Nainan [2], 
(n) the struc turc is more stable [2,5], (ii:) contact pres- 
sure distribution is better [4], and, (iv) this type of 
foundation is very suitable for the foundation on weak 
soils . 

The problem of shell foundation is useful not only in 
Civil Engineering but also in aerospace structures, solid 
propellent rocket motors, pressure vessels etc. 
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The problem of structure-foundation interaction is 
generally solved by incorporating the reaction from tho 
foundation on the struc fcure. This is done by idealising 
foundation by a simple mathematical model. Even if the 
foundation medium is complex, in a majority of cases, the 
response of tho structure at the contact surface is of prime 
interest and hence, it would be of immense help in the analy- 
sis, if the foundation can be represented by a simple 
mathematical model, without forgoing the desired accuracy. 

To accomplish this objective, many authors have proposed 
foundation models 5 a briof review of which is given below. 

1.2 BRIEF RE "VIEW OF THE POUimATION MODELS : 

More than hundred years ago in 1067 Winkler [5] 
proposed a foundation model, known by his namo as Winkler’s 
model. He assumed tho foundation model to consist of closely 
spaced independent linear springs. If such a foundation is 
subjected to a partially distributed surface loading ’ p‘ 
the springs will not bo affected beyond the loaded region 
[Fig. (1.1)]. In actual cases the foundation is observed 
to have the surface deformations beyond tho loaded region 
[Fig. (1,2)], It can be seen that Winkler’s model essentially 
suffers from a complete ’Lack of Continuity' in the supporting 
medium. Tho load -d eflec t ion equation for this case can bo 
written as, 


where lc is the spring constant and is often referred to 
as ’foundation modulus', w is the vertical deflection of 
the contact surf see and p is the load. Although this model 
leads to inconsistencies, it is hy far the simplest model. 

Even after a century this model is used for solving variety 
of problems, though their usefulness is limited. 

Other extreme is to assume the foundation medium 
as a continuous clastic solid. Using this hypothesis physi- 
cal behaviour of an actual foundation can be closely simulated. 
But the amount of mathematical complexities that creep in 
make this method unusable for practical engineering problems. 
Despite several mathematical complexities, solutions were 
presented on this basis [6,7], which, however, were limited 
to simple cases. 

The noed for bridging the gap between those two extreme 
and limiting cases and to arrive at a physically close and 
mathematically simple foundation model has been felt for 
sometime. Several authors have proposed foundation models 
which involve more than ono parame tor for the characterisation 
of the supporting medium. 
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Filoncnko-Boro dich [8] modified the Winklor's founda- 
tion by connecting the top onds of bho spring by a stretched 
clastic membrane subjected to a constant tension field 'T 1 , 
thus providing continuity in the foundation (Fig. 1,3). The 
equilibrium in the vertical direction yields the equation, 

p -- kw - T w (1.2) 

where p is the distributed normal load applied on bho surface 

of the foundation, w is normal dofloction of the surface, 

2 

V is the Laplace operator, lc and T are the bwo paromotors 
characterising the foundation. Hetenyx [9,10j achieved the 
continuity in the Winkler's foundation model by embedding, 
an elastic beam in the bwo dimonsional case, and an elastic 
plate in the three dimonsional case, wibh the stipulation 
that the hypothetical beam or plate deforms in bending 
only [Fig, (1.3)], In this case the relation of load p with 
deflection of the surface w, can be derived as, 

p = kw + D V ^ w (1.3) 

where L is the flexural rigidity of the embedded beam or 
plate, lc is the spring constant and \/^ is the biharmonic 
operator. 

Pasternak [11] suggested a foundation model by provid- 
ing for shear interaction betvree n the VJinkler' s springs as 
shown in Fig. (1,3). The shear interaction between the springs 
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has been achieved by connecting the ends of the springs to 
a beam or a plate (depending on as to whether it is a two 
dimensional or three dimensional case), consisting of 
incompressible vortical elements, which hence deform m 
transverse sheer only. The differential equation relating 
the load p with the deflection of the surface w, can bo 
derived as, 

p = lew - w (1*4 ) 

where p, is the shear modulus, lc is the spring constant and 
2 

V is the Laplace operator. It con be seen that this foun- 
dation model also consists of two parameters lc and p., and 13 
equivalent to the model proposed by Filoncnko-Borod ich[Eq. (l . 2 ) ] . 
Similar model has been proposed by Wieghardt [12]. 

Pasternak proposed another foundation model [13] 
consisting of tern layeis of springs connected by a shear 
layer In between [Pig. (1,4)]. The relation between the load 
p and the deflection of the surface w can be expressed as 

(l + “) p jr''7 2 P = kw - p-'A (1.5) 

where c and k are the spring cons Lanbs of the upper and 

* 

lower layers of springs and p, is the shear modulus of the 
shear layer. 
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In all the above models, the l/inklcr* s model has been 

modified by providing for sane interaction between tho 

spring elements and hence assuming the continuity of tho 

foundation to some degree, which is not there m Winkler's 

model. In contras b to these, starting from the clastic- 

continuum thooiy , and mbrod uc mg simplifying assumptions 

with respect to tho expected stresses aud displacements, 

some models were proposed. One such contribution was from 

Reissnor [14], who assumed that tho m-plane stresses, 

a , a and t are negligible throughout the foundation 
A y Xj 

layer. Also the hoxnzaital displacements at the upper and 
lower surfaces of the foundation layer were assumed to be 
zero. Proceeding with abo-ve assumptions and solving the 
elastic continuum equations, the equation relating the 
applied distributed surface load p and tho resulting surface 
displacement:' w, has been derived as, 


c-,w - c 0 V w = p - 


? 2 
r ^ P 


( 1 . 6 ) 


where 


E/lI and c, 


H|i/3. 1*1 and |i are tho elastic 


constants of the foundation material and II is the thickness 


of the foundation layer. It can bo seen that Eqns. (1.5) 
and (1.6) are similar. Also, for constant and linearly 


varying loads, this equation can bo seen to be mathematically 
equivalent to equations (1.2) and (1.4) hence establishing tho 



7 


similarity of tlic models. Jn this cas., neglecting the 
in-plano stresses, it can be shown that shear .stresses t 

X Zi 

end x 3 y aro cons tan L throughout the death of tho foundation 
for a given surf re o point, which is inconsistent with tho 
actual foundation performance, especially for thick founda- 
tion layers. 

Vlasov and Loontev fl5] have developed a foundation 
model starting from clastic continuum theory. Using Vlasov's 
general variational method, the load displacement relations 
can be derive d as , 

p = lew - 2tV 2 w (1.7) 

where p is the distributed surface load and w is the normal 
deflection, k and t arc the two parameters characterising tho 
foundation and con be expressed::, in terms of tho clastic 
constants of tho material and geometric properties of tho 
foundation layer. Comparing Bqn. (1.7) with 13qns . (1.2) and 
(1.4) it can bo observed that, this model is equivalent to 
tho mod ols proposed by Filonenko-Borodj eh , Pasternak and 
Wieghardt . 

ICaraoswar Ran , Das and Anandakrishnan [16] presented a 
foundation model starting from olas tic-continuum theory and 
incorporating horizontal displacements of the clastic founda- 
tion. Variational approach was used bo find the relation 
between load and foundation displaceraon bs u, v and w. 



,1 


A cl oso examination of th<~ various models ro viewed above, 
rovoals thab those models f nil short for direct application 
to practical problems, either bocauso the analysis is oumbor-» 
some or because the assumptions made for natural foundation 
media cannot bo Cully justified. More 1 over, the models 
discussed above cannot suitably be applied to shell founda- 
tion problems. 

1.3 BRIEF REVIEW OR SHELLS ON ELASTIC FOUNDATIONS; 

Some work has been reported m the literature about 
the construction of shell structures for foundation purposos 
but no work has been reported on the analysis of shells on 
elastic foundations. Shells for foundations, so far, have 
been dosigned for an assumed soil reaction. No analysis has 
boon done to find out the actual distribution of soil reac- 
tion under the shell when the shell is subjected to axisymmo- 
trie and antisymmetric typos of loads. 

Felix Candela [17] is the first to report about 
shells being used as foundation structure in Mexico, Ho 
used hyperbolic paraboloid shells. The mam reason for 
him to use shell was bo properly distribute the load on a 
weak soil. 

Shells for foundation purposes have also been used in 
India and m some other countrios. Banerjee [18] and Kaimal [19] 
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have reported works done in India using shell, rr f ounda bion 
structure. To mention a few, hypar rooting was conebruoled an 
Calcutta for building works because bhe site had soil of 
filled up nature. A water tower of 15,000 gallons capacity 
was built at Adityapur, Bihar, on a single column of 51' 
height. The base was subjected to direct load, wind load 
and earthquake moments. A hypar footing of 30 ems thickness 
covering an area of 6 m x 6 m was adopted. A dome structure 
was built for foundation for an intake well at Hatia dam. 

Hypar footings we re provided m Madras for a housing project. 

Sondhi and Patel [20j have designed hypar footings for 
a building works n n Mombase , Kenya. Conical shells have been 
used for foundations of towers and tall, chimneys at Stuttgart, 
Dresden, West Berlin, Hamburg, Hannover and Kulpenbcrg as 
reported by Ruhle [l] and Leonhardt [21], The need for shell 
foundations for tower shaped structures have been stressed 
in their papers. Using membrane analysis and considering soil 
reaction to be vortical and uniformly distributed, JCurian [22^,231 
has calculated the ultimate strength of hypar footings sub- 
jected to axisymmet ric loads. Principle of virtual work to a 
general mechanism of failure was used. 

Concrete spherical shells subjected to axisymmetric 
and anti -symmetric loadings havo been analysed • using 
membrane theory, by Sharma [ 24 J and Shanna and Mawal [3]. The 



soil reaction has been taken as umfoiuily distributed in 
axi symmetric ease and a Tune tion of Goes © m antisymmetric 
case. Using variational method, Vlasov and Leontov [ 15 ] have 
solved axisymmo t n c problem of tliin shaldow spherical shell 
on elastic foundation. They considered the normal displace- 
ment in £ 0 und a 1 1 on mo d iu m . 

Model studies have been made by many authors on shells 
on clastic foundations [ 25 , 26,27,4 ] . Nicholls and Izadi [25] 
have tested cone mod ols and hypar models on sand foundations. 
Measured strains and soil reactions arc reported. The results 
obtained experimentally were compared ;n tb the theoretical 
results obtained using membra no analysis for shells and 
assuming a uniform sand reaction. The material used for 
models were plexiglass and epoxy. 

1.4 OBJECT AND SCOPE OF THE PRESENT INVESTIGATION^ 

An attempt is made here to improve upon the foundation 
models given by Vlasov and Beontev [15] and Names war Rao , Das 
and Anandakrishnan [16] by formulating the foundation model 
in curvilinear coordinates and considering displacements in 
all the three coordinate directions. The mathematical model 
is derived using gene rails ed variational method to the elastic 
continuum. The model presented is very useful and simple for 
the analysis of any type of structure on elastic foundation. 
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In Chapter II, the equations lor the proposed founda- 
tion model has been derived m orthogonal curvilinear coordi- 
nates. Vlasov's [28] shell equations are used m formulating 
the problem. Shell foundation interaction equations are 
presented m orthogonal curvilinear coordinates, Kelnins's 
multi segment method [29] which is used to find out the 
solution of the boundary value problem, is statod in brief. 

In Chapter III, the generalised equations Cor shell 
foundation interactions are reduced for axiaymmetric thin 
shallow sphencal shells and conical shells on elastic founda- 
tions. These are then solved for normal load on shell trans- 
ferred by the column. Results are presented for static case 
and for three types of boundary conditions namely, simply 
supported, fixed and free. 

In Chapter IV, the generalised equations of Chapter II, 
arc specialised for anti-symmot nc thin shallow, spherical 
and conical shells. These are solved for radial load and 
moments transferred by a column on shell. Results are 
presented for static case and for threo types of boundary 
conditions namely simply supported, fixed and free. 

In Chapter V, conclusions and general recommendations 
for future studios have been made. 
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CHAPTER II 


FORMULATION OF THE PROBLEM .USING VARIATIONAL METHOD 

2.1 GENERAL: 

In this chapter, three dimensional equilibrium equa- 
tions of elastic body m curvilinear coordinates are converted 
to two dimensional ones by the application of the principle 
of the virtual work. In this process, the elastic body 
bepomos discrete m one direction but remains continuum m 
other two directions. This is achieved by taking the displace- 
ments as partly known in one direction and partly unknown 
m other two directions. The resulting two dimensional 
equations of the elastic body can be specialised for an 
elastic foundation on whict an elastic shell can rest. By 
the help of an appropriate free body diagram for the shell 
the equilibrium equations for the shell can be modified to 
include the elastic foundation effect. For this purpose 
Vlasov's shallow shell equations are used for the shell 
portion. 

To solve the equations of the shell on clastic founda- 
tion, shooting method, commonly known as multisegment method 
to Structural Engineers has been used and this method is 
described in brief at the end of this chapter. 
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2.2 FOUNDATION MOD 21/ ; 


2.2.1 Ba sic Equations^, for Dlastlc_ Foundation ; 


Consider a smell elmoxit of an elastic body in ortho- 
gonal coordinate system a, p and y as shown in Fig. (2.l). 
a aQ , and a are the normal stresses and cr^, cr^, and 

a etc. arc the shear stresses. let P^, Pp and P^ be the 
components of the body forces pet* unit volume on a, )3 and y 
direction^ respectively. By bhe principle of virtual work, 

( [ 

'dn- 

>, T-T \ TT !>, T 4 "ft T, 

1 


(°’/v.v ) -i- ^ (o„ c H-.E,) + <£■” H-, H 0 ) 


ba 


act "rV ' dp ^ap x± l“3' ' dy K "cty “1“2‘ 
b H_ b IL *11, 

a PP H 3 Fa" " °YY i! 2 Tv + ° Pa Ai J Fp" + °ya H 2 F7 


and 


= 0 


+ r a h 1 ii 2 '1 3 3 su drL 

^Fp ^ a pp H 3 H 1^ + 5y ^ 0 Py H l H 2^ + Fa ^ a pa H 2 ri 3 * 

biU 


(2.1a) 




^ H & I1 1 

” °YY :T 1 “5 P " G a« *3 Fp + 0 Y p il l 57 

»* «o:p ri 3 FF + P p lI lVS J 6v da = 0 


(2.1b) 


1 


(a yv H^ip t -(o Ya H 2 H ? ) + |p (« Yp 1I3HP 

F % fc n 2 

•' °aa H 2 F7 ~ G PP 1J 1 *jfT + G «Y H 2 Foe 
+ CT py H 1 Fp + P Y H 1 H 2 H 3 ^ 6w <3 - a ~ 0 




(2.1c) 
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where d-fl is the elemental volume, H-j , Hr, , are the 
coefficients of the first quadratic , 

ds 2 = I-^ 2 da 2 + II 2 2 dp 2 + H^ 2 dy 2 (2,2) 

where ds is the differential! distance and, 


u 2 '^r’ M 2 _ bjr or 1T 2 o r 

H 1 ' liDa f H 2 ~ b (3 Sp ? l S ~ S y £ y 

where r is the radius vector. 


(2.3) 


The displacement functions can he presented m a 
finite series in the form, 

1 


u(a, J3,y) 

= £ 
i=l 

(a. 

P) ^(y) 

(2.4a) 

v(a, p, y) 

m 

= ~2- v 

0=1 J 

(a, 

p)X j(y) 

(2.4b) 

w(a, p, y) 

II 

'M* 

j r 

sj 

(a? 

P)V k (y) 

(2.4c) 


k=l 


where 0 i (y),' 5t ^ ( y ) and "/^(y) are known (assumed) dimension- 
less displacement functions consistent with the constraints 
of the problem , The unknown functions of displacements 
u x ( a ? /3) s, v^(a,p) and w^ (a, p) are referred as the genera- 
lised displacements and have the unit of length. 

On substitution for displacement functions from 
Eqn. (2.4a), Eqn. (2.1a) can be written as, 



I 


[ 

V ‘ 

({ik <%' n 2 i[ 3 ) 


s 

+ ?Tp (a ap 


llg) 

0 f 


1 


an 




°ay 

n i H 2 (Hj 0 f ) - 

a 

r 

H • — — 

P 3 ft a * 






ft Ih 



■j H-i 

a yy 

H 2 vT h + v 

Li, 

b~p + ( 

y 

ya 

h 2 



+ p* h 1 h 2 ii 3 0 f ; il 3 d Y + jOayZx ^ 2 e 3 6u r dA = 0 

(2.5) 

where use of the follow on,, idonLj.Lj.ei3 have been made, 

( r, r 

j"y (°«y “lV oa i t. = [["ay 'W^V iu f dA 

-( a ny "l'Wf)’ 6u f Sa 

( 2.6o ) 

and , 

dj. = Hj dy, dA (2.6b) 

where prime denotes derivative with respec t to y. 

Since 6u.£ ^ 0, lilqn* (2.5) can bo written as, 


f, 
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Similarly Eqns . (2.1b) and (2.3c) can be written as , 

( C Sp (0 P^ ~ %y H 1 H 2 W 

a bH, ML 

+ & a (a ap H 3 H 2^ g " °yy L 1 ~&p 7 'g " °aa H 3 T p 

V H_ Mlp . 

+ a o H-. y- — X ■+• o' p II_ 's — ^ t P p H^HyJ.X ] IL dy 
yP 1 ^ g ab3dag (3 1 2 3 g J 3 1 


+ [a„ H,It 0 H, % 1 = 0 

L py 12 3 g J y 


(2.7b) 


( *, — 1 ? 2 ? * . * p* ) 


f [ - °YY Wb V'h)’ + h < V W*h 

Y % , blip 

a (a - ~ o H 0 — - Y, - or , fl PL y-yY, 

+ ^g ap31 /T h aa 2 yh pployb. 

a PI, blL 


+ ««y n 2 °py H 1 -4 y h + P Y “l^W *3 ^ 

+ [ CT yy H l H 2 H 3 ^hL = 0 <2* 7c) 


( h = 1,2,... ) 


The stress-s b ram law for isotropic Hookon material 


can be written as, 


CT ac' = + 2G e ac 


(2.8a) 


a j3p = XA + 2G e p p 


(2.8b) 


a yy - A A + 2G 


(2.8c) 


(2. 8d) 
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°&r = & 


e Py 


a = G- p 
yc: ya 


where A = = e aa + e ^ + e yv 


(2.3e) 
1 2.81’) 
(2.9) 


act’ e j3p’ ^yy 


e„„, are the normal strairsand e e - , 
"' a/ <x p py 


e^ a are the shear strois? "X and Gr are Lame's 
constants given ns, 


A = 13n/(l^)(l-2n) ) 


(2.10a) 


and 


Gr = T!/2(l+|a) 


( 2. 10b ) 


where is! is the Young's modulus and |i is the Poisson’s ratio 
of the material . 


The strain-displacement relation can be written as, 
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Stresses m terras of (generalise cl displacement functions u.^ , 


Vj and Vi c can be written as 
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On substi fcu fc ion of stresses -froi J3qnc. (2.12), Eqn (2.5) 
can be written ?s, 
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*2 £l ^ V ^ ) ^ ] H 1 ^ H 3 dY 




+ \ P y dy + r CT y y Hl H 2 !1 3 ^ h ] Y 

(h - 1,2,... ) 


= 0 


( 2 . 13 c) 
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Accuracy of the foundation model depends upon the 
number of terms taken m the series of Eqns . 2.13, and upon 
a proper choice of the functions G^Cy ) , Xj(y ) and'Y^y) 
guided by experimental results or from some theoretical 
co ns id erations . 

2. 2. 3 Choice of F unctions (2f( y ) . X- ( v ) and ~}tSxl ' 

The dimensionless functions 0(y), X(y) and "V^(y) have 
to be chosen from the physical constraints of the problem. 

For elastic foundations of finite thickness fixed on a 
rigid base, the distribution of displacements through the 
thickness can be assumed to be linearly decreasing, espe- 
cially for thin layers. In the present investigation, for 
finite elastic layers, the following functions are used, 

0(y) = X(V) = V-(r) = Sg 1 (2.14) 

where Ii is the thickness of the foundation layer. If the 
foundation layer is relatively thick, or of infinite thickness 
the choice can be made as, 

Sxnh C ( II— y ) 

0(y)=X(y) * > (y) (2.15) 

SinK C H 

where C is a coefficient depending on the elastic properties 
of the foundation and defines the rate of decrease of dis- 
placement with depth. Another possible choice for foundation 
extending to infinity can be, 
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0(y) = X (y ) = V (y) = e“° Y 


(2.16) 


where c is constant coefficient d e pend tr> g on the elastic 
properties of the foundation. Ilany such expressions can be 
selected depending upon experimen tal or theoretical 
considerations for the given nature of the foundation. 

2 , 3 SHULL E QUA! 10 MS ; 

Let a, p and y be the throe curvilinear orthogonal 
coordinates as shown m figs. (2.2). forces and moments 
acting on the shell element are shown m figs. (2.2a) and 
(2.2b) respectively. 


The equations of equilibrium car be written as. 


(I % } - n 2 |s + (as 2 ) + S x -I 


+ AB k-^Q-j_ t ABX — 0 


(2.17a) 


Sp (AN 2 ) ~ ^°1' - - 2 Za 

+ ABk 2 Q 2 + A7Y - 0 

- 0t 1 H 1 + * 2 V + j® [j| (BQ X ) + (AQ 2 )] + z 

ss (bm i 2 ) + m 2 i|I - rp (A V + M il? - “ q 2 


^A . ^ 


+ stt (BS-, ) + S. 


6) 


B 


(2.17b) 

0 (2.17c) 

0 ( 2. 17d ) 


^ < m 2 v + m i 2 % - h ( ™i) - « 2 Is - “h ■ ° <2 - i7e) 



where A = A (a, p) and B ~= ]3(oc, p) are the coefficients of 
the first quadratic I or in of the surfoc - in orthogonal 
coordinates a and (3, 

O p O O p 

ds = A da + B^d <p (2.18) 

k-^, 1^2 are "the curvature, X,T are the projections of the 
load along the pos 1 1 :i ve tangonte to the coordinate lines 
a. and /3 respectively, Z is the projection of the load 
along the exterior normal [Pig.. (?.2o)]. 

N-j_, S-^ and P G 2 are normal and the shearing forces 

acting on the area of a normal section of the body of a 
shell along the lines a - cozistant and p - constant res- 
pectively, ^12’ ^1 an( ^ ^2’ ^21’ are ^ en<3in S 

moments, twisting moments and transverse forces towards the 
positive tangent to ‘tire lines a = constant and p = constant 
respec fcively , 

The first three of the P»qns. (2.17) express the condition 
that the sum of the projection of all forces acting on the 
isolated element of the shell m three mutually orthogonal 
directions, is equal to zero. The fourth and fifth equations 
are obtained fran the condition that the moments of all 
forces acting on the isolated element of the shell about the 
two mutually perpendicular axes lying in a plane tangent to 
the middle surface are equal to zero. 
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The sixth equation of equilibrium ob turned by equating 
the sum of the moments about an axis directed along the normal 
to the middle surface to Kero , will have tb e form, 

3 1 » S 2 + k L M 12 - k 2 M 21 = 0 (2.19) 

Bqn. (2.19) can be shown to be satisfied identically with 
the use of fo rce-d cformation relation mentioned below. 

Force -d e for mo kLon relations can bo written as, 
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where strain-displacement relations are given as, 
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* W S 

A 2 B *“ U 


k l " k 2 r A £ 


rA v /^S\ B <3 ,'Sn 
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B ^ 


v 


2 r * W s _ 1 2a* w s 
AB L ^aTp A ?a~ B e>|3 


(2. 2 Oh) 

■] ( 2. 20i ) 
] ( 2 . 20 j ) 

( 2 . 21 a) 

(2.21b) 

(2.21c) 

( 2 , 21 d ) 

(2. 21 e) 


( 2.21f) 
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where u Q , v g end are the dis place laents in a, |3 and y 
directions respectively, e-^, £2 are the normal strains in 
a, j3 direc bions, 60 is the shear strain characterised by the 
change in right angle between the coordinates a = constant 
and p = constant , ^ ^2 are associated wi th bending deforma- 
tion m a and j3 directions find x represents relative twisting 
deformation . 


In case of thin shells the moment terms containing the 
factors k-^, k 2 ? k^,. k^ and the derivatives of these quantities 
and associated with the tangent displacements u g and v s have 
very little effect on the internal stresses and deformations 
of shells. The smaller the relative thickness of a shell, 
which is determined by the quantity h/J^^ where R m j_ n is the 
smallest radius of curvature and h is the thickness of the 
shell, the smaller the role of these terms. 

In the coso o f thin shells m which. h/E^^ <r l/30 all 
terms in basic equations which contain tlic quantities 
h^k-j /l2 , h 2 k 2 /l2 and partial derivatives of these quantities 
can be neglected without appreciable error. ^ 


The formulas for deformations component of a shell 
[Bqns. (2.21)] reduces to 

1 ^ A 


1 du s 


e i = Tsr + srp T s 


v - + Vf 


(2.22a) 
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1 n 2 , b 

T ST + AB ?7a u s + k 2 w s 


B 2> T 


B 

' S" 

A 

b 

B 



1 


'A 


_1_ 


(±Zls) 1 °'s 

' A 5 *5Tq ’ J>hQ 


1 *A^ W s 


(_k^!s) _ 1 *B »» 8 
1 B S|3 ' ^2-g 'ha ba 


, Tj^w 
1 r s 


1 ^ W s 


1 7>B ^ 


(2. 22b) 


( 2,22c) 


( 2. 22d) 


( 2, 22e ) 


x r ~ y x b x bi / 9 OOP \ 

AD L hadj3 “ A ^ 3a ' - 'B ®a $T J 1 ‘ ^ ' 


Bqns. (2,20) acquire the following simple form when 
small terras proportional to the quantities h^k^/12, h^k^/12 
are neglected. 


/ , S 2 \ ^ E 1 + e 2-^ 
U-^ s ) 


(2.23a) 


Eh 

27 f e 2 * !i s 

l"* |l s > 


(2.231) 


S X = s 2 = 


2( 1-t-^g ) 


(2.23c) 


" , 3 2 v &1 + F * 2^ 

12(l-p s ) 


(2. 23d) 


12(1-11^) 


+ V B § J 


(2.23e) 
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Mi 2 = M 


21 


fi h' 

... ft - 


12 ( l+(-tg ) 


(2. 23 f ) 


Q-i 


Qc 


E [ ^ um 12 ) + m 12 |f - a 4 < BM i> + M z fi ] 

(2. 23g) 

51 t ‘S’a ^ BM 12^ + M 12 " 5]3 + M 1 ^ 

(2.23h) 


liquations of equilibrium [Eqns. (2,17)] for thin shell 
can be rewritten as, 


£ (BNl) - N 2 If + ^ (as; + S + ABK = 0 (2.24a) 


(AN 2 ) - iql! +£- (BS) + S |a + ABY = 0 (2.24b) 

- (k^ -I- 1 c 2 N 2 ) + ^ [£~ (BQ X ) + (AQ 2 )]+Z = 0 (2.24c) 

^ (BM i2 ) + fl 12 |a (M 2 ) + ABQ 2 =0 ( 2, 24 d ) 

'&$ < M 12> + M 12 If " fta (DH ll + M 2i"7 " ABQ 1 = 0 (2 ‘ 24s) 

In the first two equations of (2.21), the quantities k^, 
k 2 Q 2 which result from the moments are proportional to curva- 
tures k 1? k 2 of the shell, are neglected due to their smallness. 
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By excluding from the Bqns. (2.24 ) the transverse forces 
Ql f Q 2 and expressing the stress resultants Ng? S and the 

moments M 2 ? M 12 by means of defo relations e^? e 2 J^> 

^2 and t [Bqns. (2.23)] and these deformations expressed in 
terms of the displacements u g ? v s and vr g [Bqns. (2,22)] , and 
making use of compatibility conditions? the equations of 
equilibrium can be written as [28]? 


1_ 5^o 

A da 


1 &1o 
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“ (!”l-i s ) b '~r$' + ( Ku s A 5"cT 
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(2.25a) 
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( 2,25c) 


where? 
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[j£ ( Bu J + Ip < A vJ] + 


(lc^ + k 5 ) w e ( 2.2 5d ) 



_JL_ 

2AB L £a v ^ v s' - bp 


TV [-5w7 ( Bv a ) 


(Aujj 


and , K 


r 2L /B & \ , B_ /A .2) \-j 
AD *- 1 A2g J ^ bp { B r 6 p'- 1 


^1 * ^ c 2 


(2.25e) 

(2.25f ) 
(2.25s) 


^'qns. (2.25) represent blun shells m berms of displace- 
men ts u v and w . In addition to t]ie assumptions made 
Cor thin shells, the general technics L theory oT thin shells 
is based on another assumption that the middle surface of 
the shell possesses metric of the Euclidean geometry. The 
geometric sense oC tins assumption, which is very important 
in shallow shells, is that the e:jqp res 'ions for the first 
quadratic of the surface, 

ds 2 x. A 2 da 2 + B 2 dp 2 (2.26) 


irrespective of the Gaussian curvature (A = lc-j , of this 

surface is identified with the analogous expression for tho 
first quadratic form on a plane. 

The Gauss equation, 


S'a 




d ,JL_ ^As 
dP [ B c)p J 




(2.27) 


which holds good for any surface, is D-eplpced toy the following 
simpler equation in the case of a shallow shell outlined toy 
a comparatively small part of the surfa.ee $ on the strength of 



the hypothesis it con be written, 


iL u 2. 

^ A 1d } ‘‘ 3 p 


(A- lA) _ 


(?. 23 ) 


irrespective of the value of the Gaussian curvature. This, 
means that if one superimposes on a small portion. oC the 
surface two families of lines corresponding to any system 
of or thogo nal coordinates then these lines will not differ 
m any respect from the corresponding coordinate lines on 
the plane over which the given surface rises, 

2.4 NUMERICAL iYLJ'i’nOL 

Complex at iuc tural systems neces.Mtate the use of 
approximate methods for static and dynaaiic analyses. With 
the advent of the high speed computers, the method of step- 
wise numerical integration has proved to be a powerful and 
accurate meihod for a large classes of problems. 

There are number of methods for the numerical inte- 
gration of a system of linear or dint, ly differential equations. 
Many of those methods can lie grouped under tho class of 
predictor-corrector methods or Rungc-.Cutba methods. Since 
boundary value problems are of major concern to us, the 
above methods which are suitable for initial value problems, 
should be modified to suit tho boundary value problems. 

Such methods are called shooting methods [30], These 
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methods arc capable of solving both linear and nonlinear 
bound ary value problems. Shooting mo I-pods are also called 
segmentation method and mulbisegment me tuod. In these 
methods the differential equations of m tli ordei are f ins t 
converted into a set of m first orde: differential equa- 
tions which arc two point boundary value problems. 

In case of boundary value problems since only part 
oL the boundary conditions at each bo unci 015' are prescribed,, 
a possible method of using step wise integration would be 
to start with a sot of fcrarl values for the missing initial 
conditions at the starting point. Tivm, the calculated 
terminal conditions can be compared vlUi the prescribed 
conditions at that point, and further trial can be made if 
there is no agreement between them. uowever, such a hi t~ 
and miss method is not suitable lor mixhj.no computations, 
llonco a systematic method of dote running the missing initial 
conditions is to bo used. Once the boundary .value problem 
is converted into a initial value problem, the direct 
integration approach can be used to Cmd the solution. 

While the direct integration method has a definite 
advantage, it also has serious disadvantage? that is when 
the length of the shell is large, a loss of accuracy 
invariably results . This loss of accuracy does not result 
from accumulative errors in integration but it is caused by 
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the subtractioi oj almost e qu d n umbers in the process of 
(-CD en ' in. • lioi of B . o unlm own bounds r values. Tha t in, 
for every jet of geometric and materiel parameters of the 
shell there is a eritical length bevond vtIiicIi the solution 
looses al i accuracy. The advantage oj. finite difference 
approach over diiect integration is that it can avoid such 
Josr. of a ecu " . It is concluded f ?9J that if the 

■ solutions of Lh“ .ivstei o ul* obric actuations, which results 
from the finite ci l Pf erono c equation i. li obtained by means of 
(reus -:ien elo miiio ti on , then no loss of ax curacy is experienced 
if the length of the shell is increased. 

The method used hero is giver by Jtalnins [29]. This 
method lias definite advantage over the finite difference 
approach. The main advantages ere (a) it can be applied 
convonientl ' r to a largo system of first oi‘der differential 
equations and (b) it permits selection of an optimum step 
sj 20 of Integration at oneb step according to the desired 
accuracy of the solution. 

The method used hero can be divided into two parts 
(a) direct integration of (m + l) initial value problems 
over preselected segments of the total interval, and (b) the 
use of G-aussian elemmation affor the solution of resulting 
system of matriA equations. 
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The initial volue problems are defined over segments 
of the total m towel, the length of which are within the 
range of applicability of the direct integration approach. 
After the initial value problems are integrated over these 
segments, con tin.n ty conditions on all variables are 
written at the end points of the segments and they constitute 
a simultaneous system oi‘ linear matri’ equations. This 
system oi mat nr equations is then solved directly by 
means of O-aussian elimination. The result is that direct 
integration approach is applied and there is no loss of 
accuracy because toe lengths of the segments are selected 
in such a way thane the solutions of the initial value 
problems are kept sufficiently small. 

2,4.1 lie d ue ti on of Boundar y Vslne P rob lem to Initial 
Value Pr ob 1 em : 

Two point boundary value problem governed by differen- 
tial equation, 

= A(x) y ( x) + B(x) (2.29) 

dx 

is reduced to a series of initial value problems. In (2.29) 
y(x) is an (m,l) matrix which represents m unknown func tions 5 
x is the independent variable,* A(x) denotes the (m, m) 
coefficient matrix and JB(x) is the (m, l) matrix of non- 
homo^neous terns. The objec t is to determine y(x) m the 
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interval a ^ x < b subjected to m boundary conditions in 
berms of linear combination oi y(a) and y(b) in the form, 

I'a y(a) + P]) y(b) = 5 (2,30) 

where F , ar (m, m) matrices and G is an (m, 1) matrix 
which are known from the statement of the boundary conditions 
of the problem . 

Jjeb the complete solution of (2.29) be written as, 

y(x) = f(x).C + Z(x) (2,31) 

where the (m, 1) iia tnx 0 represents m arbitrary constants 
and Y(x) is an (m, m) ond Z(x) an (m, 1) matrix which are 
defined as the homogeneous and particular solutions of (2,29) 
in the form, 

= A (x) Y ( x) ( 2 . 32 a) 

and £§^1 - A(x).Z(x) + B(x) (2.32b) 

The initial conditions for determining Y(x) and Z(x) are 

Y(a) = I (2.33a) 

and Z(a) = 0 (2.33b) 

where I is a unit matrix . 
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evaluation or (2,.;1) at x = a leads at once , m view 
ol Eqns . ( 2.33a) , ( 2.33b) to G = y(e) and then (2.31) at 
x = b can be written as, 

y(b) = 1(b) y( a) + 2(b) (2.34) 

Together with (2. r '0), equation (2,34) constitutes a 
system of 2m linen r al gebrare equations from which the 2m 
unknowns y(a) and y(b) are determined. Once y(a) is known, 
the solution at -my value of x is obtained from (2,31) 
provided Y(x) and Z(x) at that particular x are stored. This 
completes the reduction of o two point boundary value problem 
defined by (2 .29) to (m +1) initial value problems given by 
Eqns. (2.3 2) and (2.33). 

2.4.2 Mu lti segme nt M ethod of Integra tion : 

A differential equation system of even order, say m, 
with m/2 known boundary conditions at either boundary will 
be considered. however? the method is equally applicable 
to a more general case. Let the shell be devidecl into 
M-segments, denoted by Si? i = l,2,,.,il, of arbitrary 
length. Let tho coordinates of the ends of the segments 
be x = x^, where the left hand edge oc the shell is x^ 
and right hand edge is at x = X M+1* Tiien solution of 
equation (2.29) following the analogy of equation (2,30) can 
be written as , 
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y(x) ■= Y t (x) y(x x ) + ^(x) (2.35) 

l = 1,?,... M 

where Y^x) denotes the fundamental matrix solution for the 
segment and Z is the corresponding particular integral. 

Y and Z^ are defined numerically at finite number of points. 

The intersegment continuity condition can be written as, 

y (x i+i> = y i (x 1 +i ) y(x i ) + 2 i ( \+i ) (2 * 56) 

l = 1,2 , , . „ M 

The continuity condition between two successive segments are 
written in a matrix, form for all the segments. 

.■porpr example, it can be written xthat , 

y 2 = Y i Y i + z l 
y 3 = Y 2 y 2 + Z 2 

= Y 2< y i y l + Z l> + Z 2 
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Once is known in terms of > tlie resulting algebraic 
equations can be solved using Gauss g lamination, The matrix 
can be written as. 
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( 2 , 30 ) 



Y 

( mxm ) 

= IVVl • • • • y 2 T h 

(2.39a) 

and 

Q (l) 

= Tl " 1 Zj. + Y^ 1 Yg " 1 Z 2 + ... 



(mxl) 

• •• + Y 1 1 Y 2 1 ••• Y fl 1 Z, [( 

(2.39b) 

and, 

q (2) 

'Q 

is the R.H.S. of the boundary conditions 

at 


X 1 and X M+1- 


This gives all the functional values at and 
x i v i+l* Rungc-Kutta method of integration, values of 

functions for all x in a <L x b can be obtained. 
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CHAPTER III 


AX1SYMMSTRIG PROBLEMS 


3.1 GENERAL; 

Equations for foundation inodel and sheila developed 
in Chapter I?,. arc specialised for thin shallow spherical 
shells on elastic foundations and than shallow com cal 
shells on elastic foundations, applicable to axisymmetnc 
problems . 

Shallow spherical and conical shells subjected to 
normal load transferred by the column have been solved for 
fixed, simply supported and free typos of boundary conditions 
'the foundation thickness has been considered as finite m 
all the examples that have been solved. Foundation of 
infinite thickness, layerod elastic foundation could also 
be considered without any difficulty , Results have been 
presented in non -dimensional form so that they can bo used 
for design purposes. 

In case of spherical shells two parameters, which 
accounts for thinness of the shell (h/R) and shallowness 
of tho shell (R/a) have been considered. Results have been 
presented for different values of those parameters. In case 
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of conical shells results arc prescribed Cor different values 
of 9, where 0 Is the angle made by generators and axis of the 
shell. 0 accounts for the shallowness of the conical shell. 

3.2 EQUATIONS POE THIN SHALLOW 3HE1IS POR AAI SYMMETRIC 
PROBLEMS s 

3.2.1 Spherical Shell 1 ; 

In case of shallow shell, polar coordinate system 
can be taken as, 

' ' ds 2 = dr 2 + r 2 dp 2 (3.1) 

Comparing Eqn. (3.1) with the first quadratic equation 
of surface m curvilinear coordinates, Eqn. (2.24), it can be 
seen that, 

A = 1, B = r, da = dr and d|3 = dp (3*2) 

Por spherical shells both the curvatures are equal , [Fig. ( 3 .la) J 



This gives Gaussian curvature as 

K = k, k ? = ~~ (3.4) 

ir 

4n axisymmetric problems, the displacement functions are 

independent of p and the displacement function v a m the 

s 

circumferential direction is zero. 

% = U s (r) 

w s = w s (r) 


(3.5a) 

(3.5b) 
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Substituting Bqns. (3.2), (3.3), (3.4) and (3.5) m the shell 
equations in curvilinear eoordmo leu [bqn . ( 2 . 25 )], the result- 
ing expressions are, 
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(3.6b) 


Equations (3.6a) and (3.6b) represent thin shallow spherical 
she’ll applicable only to axisymme brie problems. 

3.2.2 Conica l. S hell s ; 

Polar coordinate system can also be taken in case of 
conical shells, hence Eqns. (3.1) and (3.2) arc also valid in 
case of shallow conical shells. For conical shells [Fig. 3.1b] 
radius of curvature R^ is infinite and 


R. 


- r/Cos Q 


(3.7) 


als o K 


0 


(3.8) 
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I p or axisymmetno problems , 



u s 

- u 0 ( r ) 

>~J 

(3.9a) 


w s 

= w s (r) 

(3.9b) 

and 

v s 

=■ 0 

(3.9c) 


Substituting Eqiis. (*5.2), (3.7), (3,8) and (3.9) m the 
shell Tiqns . (2,25), the resulting expressions are, 


+ 1 % U S Cos © , . ° 0S ® dW 3 

7 + r dr - J - ’’c + r TF 


dr 




ki 

ir h ' yL 
s 


(3.10a) 


.2 w 9 c] J w ( 1-Cos 2 0) d 2 w 

" 12 l;-T" + _ F "T~ * p — 

dh r dr 5 r dr 


. (1-3 Cos 2 0) dw s 4 Cos 2 0 ir t Cos 2 6 „ 

dr r 4 s J r 2 s 


l l a 0 oM 6 _ C os 9 

r dr 


(1 - ) 


u 


r 


B, h 


Z (3.10b) 


ISqns . (3.10) represent thin shallow conical shell applicable 
to axisymmetno problems only. 

3.3 EQUATIONS FOR ELASTIC FOUNDATIONS APPLICABLE TO 
AXISYMMETRIC PROBLEMS: 

Equabions for elastic founds bion are developed m 
Chapter II. These foundation equations are particularised to 
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suit for shallow spherical shell and shallow conical shell. 
'I'he first quadratic equation is taken as, 


2 2 2 ? 2 
ds * dr* + r dp + dy 


(3.11) 


Comparing Eqn. (3.1l) with the first quadratic equation in 
curvilinear coordinates fEqn. ( 2, 2 )] , it can be seen that, 


da = dr, dp = dp, dy = dy, = 1, H 2 = r 

and H 3 = 1 (3.12) 

Substituting Idqns . (3.12) m foundation equations m curvi> 
linear coordinates Eqns. (2.13), the resulting expressions are, 


2 

1 d u. 1 a „ du 

5~ a i .l 5^ i 

i=l 1X dr* i=l r ar 


1 ^ 1 

2. ~y ”i - Z- b if u i 
i=l r i =1 


n dw, n dw,. _ 

+ d kf '“dr " ^ c kf dr" = " *rf 


(f = 1,2,. 


n 

I 

lc=l 


,2 

d w 


’ kh dr" 


k i 
V + 


1 du 

- Z 
1=1 


ih dr 


y- fkh 
kTl r 

. 

dr 

n 

z 

k=l 

e kh 

1 

du. 

1 

c , i 

+ .1 °xh 
1=1 

-r- + 

dr 

z 

1=1 

ih 

‘V 


w. 


u. 


1 d.. 

y — u. = - p 

4 r l 
i=l 


where 


yh 

(h = 1,2, ... ) 


a if 


= 


•H E 


'o (1-|A 0 ) 


21 <*i <*f d Y 


(3.13a) 


(3.13b) 


(3.14a) 
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b .. 
.it 


/Li |j] 

— 0' cy 

2(1 + ji Q ) 1 J 


'kf 


'ih 


lrJ? 


-f 


H 

o 

H 

o 

II 

o 


E 


V, . 0’ dy 


2(1 + H 0 ) 


2(1 + |i 0 ) 




i ' h 


C / 


(T~77T < -r 


/IT 


lh 




o 


(1 - n 0 ‘) 


0 2 s^'h ^ d t' 


kh 


II 

"'0 


J3 


o 


2(1 + H 0 ) 


VkK t ’ Y 


/H 


’kh 


(l - ti 0 e ) 


" YlVl cly 


2s ' krh 


p 

J rX 


p r y) 0 g - d r + 0 r ]; J 


/II 


■Yh =1 P Y v) ^ dy + [a YY V h ]H 


and , 


E 


E 


(1 - (i") 




Ji. 


o 


1 - |i 


(7.14b) 

(3.14c ) 

(7.14d) 

(3.14e) 

(7.141) 

(3 .14g ) 

(3.14h) 

(3.141) 

(3.14 j) 

(3.15a) 

(3.15b) 








J C- 

cl-. 


C H 

E 

0 

01 2 C 'Y 

( 3 . 18 s ) 

l ll 


)o 

(1 - u 0 2 ) 

"K 

= { 


E 

0 

dY 

( 3 . 18 b) 

D ll 


\ 

) 0 

2(1 + hy) 


1 

( H 

1 

li 

0 

\ #1 ( 'Y 

( 3 . 18 c) 

G ll 


)o 

2U -1- |i 0 ) 

d 

_ { 

Ai 

E 11 

0 1 0 

V I h ^ 

( 3 . 18 d ) 

C 11 


)o 

(1 - n Q 2 ) 




E o 

Vi dy 

( 3 . 18 e) 

r ll 


'o 

2(1 + H 0 ) 

S 11 

_ ( 

rH 

B o 

4/ ' 2 , 

h cly 

( 3 . 18 f ) 

~ 1 

1 0 

(1 - H 0 2 > 

*'rl 

. 1 

rH 

i 

1 0 

f r (r, y) 

0 X (y) dy + [o rY (y ) 3 q 

( 3 . 18 g ) 

r yl 


'II 

' 0 

f y (r, y) 

y (-,) dy + [ c, Y Y ’V 1 (l)]o 

( 3 . 18 h ) 


In the case of externally applied loads on the surface of the 
body (along y =0), first part of the right hand side of 
Eqns . ( 3 . 10 g), ( 3 . 18 h) have bo be token as Sfcieltjes Integrals ? 
as given below ; 

•H 

i P (r, y) 0 , (y) dy = P (r) 0 ,( 0 ) 

Jo r 1 r 1 

+ | E r (r, y) (y) dy (3.19a) 
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^ P Y ( r > y)W-l (y) dy = ? Y (r)V 1 (0) + 

+ £ P y (P ? y)V 1 (y) dy (3.19b) 

where P r (r) and P (r) are the externally applied surface 

loads acting along r and y directions-; respectively . The body 

force components P (r, y) and P (r, y) are neglected. 

i Y 

Therefore , 

(H 

) p r ( r > Y) ^ (y) dy = i’ r (r) ^ (0) (3.20a) 

^ p y (r, y) 1 !^ (y) dy - ^(r^ (0) (3.20b) 

3.4 SHELL .FOUNDATION INTERACTION EQUATIONS 

3.4.1 Thin Shallow Spherical Shell o n Elastic Foundation : 

Consider a thin shallow spherical shell on. an elastic 
foundation (l |l ig.3.2a) carrying axisymmetric loads. It is 
assumed that the friction between the shell and the founda- 
tion as negligibly small. Hence the shell transfers only 
normal reaction on to the supporting medium. Thus the 
differential Eqns. (3.6) for thin shallow spherical shell 
can be written as, 
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? ?U B + i 1% _ ^5 (1 ' + (1 + “ S ) flW S 

. r dr 2 ^ 


dr 


ll 

E h 
s 


H dR" 


(3.21a) 


2 2 

. d !•/ 0 d w a n dw 

n r a + s 1 . s. _2 s 1 s 

~ 12 L . 4 r,3 2.2 h „2 . 2 h 3 <3r 

dr r dr R dr r 


.2 d^vi 
hr S 


2 fl \ 


dr 1 


r\ dW — N « 

, _2 S i 

+ o J - --- 


rR 2 dr 

U + n a ) 

3?r 


2(1 + iO (1 + |A fl ) du 

vr . - 


R 


R 


a ' _s 
dr 


u. 


(i-n 2 ) 

■ • srfr- t z - V r)J (5- 


Where P^(r) is the normal reaction due to elastic foundation 
on the shell. 


As the shell is assumed to transfer only normal forces 
to the founds bion, and assuming the compatibility of normal 
deT-lcc Lions of bho shell and foundation surface, P^(r) can 
be eliminated from Eqn. (3.21b) using Eqns. (3.20b), (3.18h) 
and (3.17b). Considering that P r (r) equals to aero in 
Jiqn. (3.17a), tlic governing equations for the displacements 
of the shell foundation system can be written as [choosing 


^(o) = Y 1 (o) = 1] 

a, n du. 


2 

a^u 


a 


11 


dr 


* + -ii ra 

5 r dh 


£L 


li 


u., - b. u. + d 


dW-, 


2 1 U 11 1 11 dr 


- c 


dWj 

11 dr" 


If 


- ~K Y 


(3.22a) 


21b) 



n 


,2 

d u 


dr ‘ 


1 du s 


u s ( 1-lip.) (l+|a Q ) dw n 


T ~dF ~ "? + 

r 


R 


T 


TT 


d r 


(l~P- e ) 


(3.22b) 


2 

d w-, r-,., dw-, du n du-. 

- n -l vj- i* n — S - W * d n n — — + Ch-i j 

If c1 2 r dr 11 1 11 dr 11 dr 


u u 

+ c n — - d n . — 
11 r 11 r 


H s h 2(1 +|A ) 


(l“ho ) 


2 X if 


o 


h 2 ( d ^ v l 


12 


+ 


.dr 


_2_ rl y _l 

r 7r T ' 


,2 ,2 

1 W 1 . £. ° V 1 

22 0 2 . 2 

r dr I? dr 


* _i n * j. fii 

r 3 dr rR dr 


)■ 


2( 1+p ) ( l+M'o ) du 

i » T * , b ^ 

'■ ~ W 1 + “~R "d7 


(l+lO n 2(l+|i n ) ( 'H 

+ R? = - “ E [Z+ f fyy W Y } ! ] 

o ( ' 1 i o 


where a-j^, b-^ etc, are given by JjJqns. (5.18). 

When free boundary condition 1'or the shell is consi- 
dered, the flank (portion of foundation beyond shell founda- 
tion) has to be taken care of [Pig. (5.2a)]. If u^, w^, are 
the displacement functions for the flank portion, the govern- 
ing differential eciuations for a free shallow spherical shell 
on an elastic foundation would be hqns . (5.22) taken together 

with the following equations which accounts for the flank. 
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a ll 

:7 + 
2 + 

d v d 

a u 

r 

dr" 

a n 

r 2 U f 


- 

G 11 

dv f 

dr 

= - [(cr 

r ll 

d 2 v 

7? + 

dr 

r il 

r 

in 

dr 

3 11 w f 


+ 

C 11 

.. 

r 

U f 

d ni -i « 
11 r 

11* 

h ll eic 

. are given by '!< 


dir 




H 

Q 


(3.23a) 




du 


+ c. 


ill 


Y Y ' 


3. 4 o.l,l Choice of Dimensionless Displacement Functions; 


In alt the numerical examples that are considered, 
it L k - assumed that the elastic foundation consists of single 
elastic material and is of finite thickness H. The dimension- 
less functions 0^(y) end y^(y) have to be chosen from the 
physical constraints oC the problem. In the present investi- 
gation, for finite clastic layers, the following functions 
arc used, 

0l(y) = ^(y) = -p- (3.24) 


Substituting 0 1 (y) andV 1 (v) from Eqn. (3.24) m Eqns.(3.18) 
for a 11? etc. and integrating, the resulting expressions 
are , 


Oi 


11 



H 

3 


(3.25a) 


j 7 


b 


11 


C 11 


d 


11 


r 


11 


3 


11 


0 

i 

2 0- + M 0 ) 

'u' 

D 

0 

i 

2(1 + H 0 ) 

2 

13 u 
o ' 0 

1 

1™ ' 

i - 

I 1 

! ij. 

2 

E 

0 

H 

20 + l-o > 

3 

B 

q_ 

1 

& - i‘ 0 z ) 

H 


(3.25b) 


(3.25c) 


(3.25d) 


(3.25e) 


(3.25f ) 


In order to non -dimens lonalise the XJqns. (3.22) and (3.23) 
the following non-dimennionalisec! functions are assumed, 
u a - u o /li j = u^/h, w-^ = T/-| /h , Uj. = w^/h, 

Wp = Wj./h, r = r/a (3.26) 

where a is the radius of the shell and h is the thickness 
o C the shell, 

St roes ou o , a arc zero a l the surface of the 
ry 1 yy 

foundation and flank and at y = H, 0^( y ) = Vy(y) =0 
[ Lfiqn . (3.2/))]. Prom above it can bo seen that, 

- KyVl<r< = [ K y >f 


[(q Y V\nn)]‘ ; 


0 


(3.27) 
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If the shell is subjected to only normal load trans- 
ferred by the column, coming as ring load, then the components 
of the forces X, Z are all equal to zero. 


Taking care of the facts stated above and substituting 
for a-^» bqq etc. from Eqns . (3.25) and also using the dimen- 
sionless functions given by Bqns. (3.26), the differential 
Bqns. (3.22) and (3.23) can be written in dimensionless foira 
as given below. 


Differential Bqns. (3.28a), (3.28b) are for the flank 
portion and are used only in case of free shell along with 
other equations, 





r dr 




,a N 2 

( H ) u. 


1 


(3.28b) 



(3.28c) 
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= 

dr 


r dr 


(1 + !'J (|) 


(1 - |X„) (|) 2 » 


(3.28d) 


s (1 + |i Q ) n K ^ a dr 2 

+ 1 (IJ) L 2 ,cu - 

r dr (l-.i ) H 1 


- rrr-jrj 




r dr- 


du l + 1 

3 u i 
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, 1 U 1 
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^ — 4 

or 

h- 2 ~ " ’ 
r 

i d2 ' w i 

o 

fS.\2 a \ 

? 2 dr 2 

— £ 

C\J 


( l-p 0 ) r 


Un 

~ ] 


2 /ax2 1 dw l 1 d¥ i { . /R\ 2 

2( k ) - =3- ™- - 24 (1 + p. ) ( E ) 

r dr r dr 


X (|) 4 ^ - 12 (1 + | 1 0 ) (f ) 2 (f ) 3 

dr 


12 U + !*„ ) < E > 2 (|> 3 Uj 


(3.28e) 


3. 4. 1.2 Forces and Momenls m Shell: 


Forces and moments in dimensionless form in spherical 
shells for axis yvmnet nc loading can be ''rritten as, 

iT. n du „ u 


n , R L 9 

(1 " Hg ) 


V, n dU t> u 

(|) [(f) ~ + ?! + H s (f) T? 

K 8 dr 1 s a r 


+ >* a »! ] 


(3.29a) 
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(3.29b) 


(3.2Qc) 


] (3.29d) 




du., 


i au i ui n 

+ 2 — + 2 — - 

dr r 


u- 


0 ' D 


where , 


D = 


R s h' 


12(1 - \i/) 


(3.29f) 


3.4.2 Thin. Shallow Conical Shelia on. E las tic Foundations i 

Proceeding on the similar lines, as m case of spherical 
shells, the governing differential equations in dimensionless 
form for the conical shell on elastic foundation can be 
written as given below [See Pig. (3.2b)]# 



Differential equations for flank portion are 
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Differential equations for shell~foun.de tion system are 
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(3.31c) 


where Uj., w^ eLc. are defined by Eqn, (5.26). 

Differential Equations (3' 31) sire solved for fixed 
and simply supported boundary conditions. Differential 
equation (3*31) and (3.30) together are solved for free 
boundary condition- 

Forces and momenta in dime ns ion less fo rm in conical 
she 11s for a xis y mmet nc 1 0 ad 1 ng can be w rit t en as 
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+ O — + P ™ 

dr dr 


2 

(i - ^ 0 ) 
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(3.32e) 


where , 


D 



12 (1 - !i s 2 ) 


(3.32f) 


3 . 3 BOUNDARY CO NDIT IONS : 

Three important types of boundary conditions namely, 
simply supported (b.S.),, fixed and free, have been considered 
for both conical and spherical shell, 

3.5.1 Boundary Co nditions at t he Outer Bou ndary o f Shell, s 
3. 5*1.1 Simply Supported? following are the boundary conditions. 


(a) u >1 = 0 (3.33a) 

(b) u - 0 (3.33b) 

(c) w x = o (3*33c) 

(d) M-j - 0 (3.33d) 


In dimensionless form, the above mentioned boundary 
conditions for both spherical and conical shells can be 
written as, 
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(a) 

Tl l * 

0 

(3.34a) 

(h) 

% = 

0 

(3.34b) 

(c) 


0 

(3.34c) 

in 

c 

d 

+ 

dr 

11 dw T 
— b - = 0 

p dT 

(3.34a) 


3. 5. 1.2 fixed or Built-in: 

following would bo the boundary conditions for fixed 

case , 


(a) 

U 1 = 

0 

(3.35a) 

(b) 

u s = 

0 

(3.35b) 

(c ) 

W 1 = 

dv. 

0 

(3.35c) 

(d) 

1 

dr 

0 

(3.35a ) 


In dimen si on less foim, the above mentioned boundary 
conditions for both conical and spherical shells can be 
written as. 


( a) 

u l “ 

0 

(3.36a) 

(b) 

-B = 

0 

(3.36b) 

(c) 

»1 = 

0 

(3.36c) 

(d) 

dw^ 

0 

(3.36a) 


dr 
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3.5 . 1.3 Free; 

(1) In free type of boundary conditions the flank 
portion of the foundation needs to be considered. The 
boundary conditions at the far end of flank can be written 


as 

(a) Wj. ■= 0 (3.37a) 

( b ) Uj. = 0 ( 3 . 3 7b ) 

(2) The boundary conditions at the free end of shell 
can be written as, 

(a) Uj-> = u^ (3.38a) 

(b) Wj. = w-^ (3.38b) 

(c) Transverse shear m flank 

= Transverse shear in foundation 

-h Transverse shear 111 shell (3.38c) 

(d) Normal force m Elank 

- Normal force in foundation (3.38d) 

(o) Normal, force m shell =0 (3.38e) 

(f) = 0 (3.38f) 


Above mentioned boundary conditions in dimensionless 
form for spherical shell can be written as, 

(l) at the far end of the flank 

(a) -- 0 

(b) w f = 


0 


(3.39a) 

(3.39b) 
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(2) at the free end of the shell 


( a ) \x.£ = u-^ 

( h ) Wj 
Tj r 

( c ) ]jr 


= W, 


6 ( l-jig 2 ) 


(3.40a) 

(3.40b) 
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(e) 

S_ y 

d r 

t|) 

w a + >- l s + (f) 
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] 

= 0 (3. 

. 40e) 

U) 

d 2 w, 
’=? + 

ta 

dv?i 
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. 40f ) 


dr 

r 

dr 






Boundary conditions in free case for conical shell an 
dimensionless form can be written as 

(1) at the Car end of the flanh 

(a) u f = 0 (3.41a) 

(b) w f = 0 (3.41b) 

(2) at the free end of the shell 

(a) Up = u-j_ (3.42a) 
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+ ::r ~ b n ™- + | (u) w. 


U., 


dr- dr 

- 1 i- l o ( i> ®f = 0 


du u Wn Cos 9 

(e) -=§•+ n [ ^ + -U_ — - J 

ar k r r 

2- 

d w, |J. , dw-, 

( f ) _i + —L. = o 

dr^ r dr 


2 |J, o V H' 


= 0 


(3.42c) 


(3.42d) 


( 3 . 4 2e ) 


( 3 . 4 2f ) 


3.5.2 Bo undary C o nditions at the Column End s 

Boundary conditions at /'the column end of the shell 
can be written as, 

(a) Transverse shear in shell 

-i- Transverse shear in foundation 

= Load/ uni t le n ;;th (3.43a) 

(b) \x x « 0 (3.43b) 


(c) u 0 =• 0 


(3.43c) 



r~ 

_ 


dvr . 

(d) 'dr " 0 (3.43d) 

In non-dimensional form, the above mentioned boundary 
conditions can be written as ; 


(l) Spherical shell 



s 


where 1 ? is the load transferred by the column on shell 
lound at ion , 


w 

u i 

= 0 



(3.44b) 

(c) 


=; 0 



(3.44c) 

(d) 

dWj 

dr 

-- 0 



(3.44d) 

(2) 

Conical she ll 




(o) 

14 

o 

) a 3 

U-j 

Ci-( 1 0 > 

t 1 ^1 

n tt dw" 

+ 1 (f ) ] 

J cl dr 




f fli A. 

d r 

d 2 ^ 

dr 2 

1 d v i 
___ ] 
r dr 




It 

- t(|H; 

3 hL 
s 

(3.45a) 
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(b) u^_ - 0 


(3.45b) 


(c) u = 0 (3.45c) 

dw-, 

(d) ~ = 0 (3. 45<3) 

dr 


In both spherical shell and conical shell? 

M Jp 

r(i)( _)J has taken as the load parameter. 

L b ( ^ 

s 

3. fa NUMERIC AL CALCULATlOikj A1TU GONCLiJSIOhS : 

In all the numeric;.! examples that have been solved, 
following parameter values have been kept constants 

U - 2.0 

! - °- 2 
H S = 0.25 

H 0 =■ 0.L7 

Jn case o(' spherical shell, (~) - 4.0 has been kept 

constant while (f<) and (^r ) have been varied. Also 

^ O jo $3 

(I 1 ) = 0.0333 has been kept constant while (") and (g~ ) have 

K ' o 

been varied. 

Similarly m case of conical shells (j|) 

E 

constant while Q and g* take different values. 


14.0 is a 
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Results have been presented in. plots drawn between 
•” and functional values at unit load parameter, load 
na-' ano ter has been taken as, 


p 

=3 

(|) (- 1 — ) 

V 

(3.46) 

functions’ have been redefined as follows s 


r*- 

U I 

- 

■ - ■ 

V J/ 

(3.47a) 

-“v 

s 

=: 

u /V 
rr 

(3.47b) 

W 1 

=■ 

w 1 /P 

(3.47c) 

"\ 

N 1 

=- 

*\- 

•r j 

t'll/P 

(3.47d) 

Wo 


V P ‘ 

(3.47e) 

ll'l 


V* 

(3.47f) 

M 

u 2 

- 

u 2 /v 

( 3 . 47g) 

r i 

q 

- 

q/£ 

( 3 • 4 7h ) 

Diinons i on Tear 

<*i 4 ^ / 

; d i epic icemen Ls and forces u-j , u g , 

*1- 

n 2> M 1? 


lip and Cor shallow spherical and conical shells on 
elastic foundations have been presented for simply supported, 


fixed and free tyjjes of boundary conditions. 

3,6.1 Simply 3uppq rt_ed_£ind Fixed Bound any. J^ojidition J. : 

(a) Spherical Shells 

From Figs. (3.3a), (3.3b), (3.3c), (3*4a)> (3.4b) and 
(3,^c) it can be seen that u 1? u g and w-^ increases with 


n 


decrease in (~ ) , also the values of junctions decrease with 

Pi 

decrease m (7^-) i.e. when foundation is stiffer. One important 
j! 'o 

observation that emi be made from Figs. (3.3a), ( 3 . 3 b), 

(3.4 a) and (3.4b) is thab the radial displacements in shell 
are much larger than those in foundation. This happens 
because the shell is free bo slip over the foundation as the 
friction between them is negligible. Ib can also be observed 
that the maximum values of and u Q occur at samo P for a 
particular parameter. It can also be observed that the 
values oJ u-^? u a and for fixed boundary conditions are 
less than those for simply supported boundary conditions. 

Responses Nj , Np for simply supported spherical shell 
ore shown in Pigs. (3.3d) and (3.3®) where it can be seen bhat 
is negative and 1S positive for all values of P and 
they increase with increase of (^) value. Pigs. (3.3f), 

(3 .3*5), (3-3h), (3.4d), (3.4®) and (3. /If) show the variation 
of’ M-j , Mjj and q along P. It can be observed that thoir 
values Jncreo.se wibh increase in (^) ratio. It can be 
no bod that v a, nab ion. of q is not constant; or near that, q has 
maximum value noor the column end and decreases away from it. 

R 

The responses u^> u s mcreoso with decrease 111 (~) 
as shown m Pigs. (3.5a), (3.5b) where as w^ increases with 
increase in (§) [Pig. (3.5c)]. u s is greater than because 
of the slip between shell and foundation. ^ increases as (-) 
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increases [Fig. (3.13d)] ’out reverse is the case for T'Tp 
[Fig. (3.5e)]. As the shell goes deep (as (~) decreases ), 
the stretching of the shell is more, and since the boundaries 
are not allowed to move, more negative Mp ano positive q 

develops near the shell boundary [ Fiqs . ( 3- 5f ) » (3 . 5g) , (3. 5 h) , (3. 6 a) , 
(3.6b), ( 3 . 6 c)]. Also the absolute valr es mf Mp and q 

decrease with increase m (r), except near the column end 

ci 

where values for all ( — ) are almost seme. It can be observed 

ci 

that tho values of Mp and q are greater for a greater 

E 

(■—) ratio. 

JJj 

o 

(b) Conical Shells; 

In case of simply supported and fixed conical shells, 

Wj decreases w± bh decrease m the value of 0 , which means 

that normal deflections are decreasing as the shell goes 

deeper [fig. ( 3 . 7 a) and ( 3 . 0 a)]. When 0 = ^ 1.57, bhe 

conical shell becomes a. flat circular plate and w at 

0 s 1,57 is as largo as 10 tunes of that at 0 = 1.31, as 

shown in Fig. (3.7s). It can also be observed that w is 

higher for higher 4 - ratio which means that w is having 

■ L 'o 

greater value when the foundation is loss stiff. 

Responses M-, , Mp and q are having much greater value 
for plate than m shells and they decrease with decrease in 
0 value, as can be seen from Pigs. (3.7b), ( 3 . 7 c), ( 3 - 7 d), h;c 
(5,8b), ( 3 . Sc) and (3.8d). 
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3.6,2 F ree Boundary Conditions : 

In case of free boundary conditions, the flank (portion 
of foundation beyond the shell) has been considered. Results 
have been presented for flank widths equal to 1.6 times 
and 2.4 times the outer radius of the shell. From Figs, (3.9 ) , 
(3.10) and (3.11) it can be observed that bhe responses 
coincide for the two flank widths con aid ere d. From this it 
can be concluded that the flank width equal to 1.6 times 
bhe outer radius of the shell ns sufficient for solving any 
free boundary condition problem. 

( a ) Sph e ric al Shel 1 ; 

From J?igp. (3.9a) and (3.9b) ifc con be seen that u g 
and 'w^ increase with the increase m (g) ratio. In Fig,. (3.9b) 
the dotted lines and firm lines represent normal displace- 
ment when width of the Clank has been considered as 2,4 times 
and 1.6 times the outer radius of the shell respectively. It 
can be observed that the normal displacements in shell 
portion hove the same value for both chc flank widths 
considered. A kink mthe curves at the boundary of flank 
and shell can be observed which is duo to the change in 
slope. is slightly negative near the free boundary but 

becomes positive as it approaches the middle of the shell and 
its value increases with increase in (jj) [Fig. (3.9c)]. 

Similar is the cose for R 2 except thab it is slightly positive 
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near the freo boundary and becomes no otive as it approaches 
column end [Fry. (3.9d)]. Figs. (3.9 o) and (3. 91) show M-j_, 

Mg lor diffe3 enfc values of (g), and they arc maximum at the 
least value of (^) . From Fig. (3.9g) ? the maximum value of 
q is same for all the values of (fe) but the curve flattens 
as (^) increases. From Fig. (3.10a), it can bo seen that 
w, mcroascs with increase m (— ). From Fig. (3.10b) it can be 

-L 3 

ru' 

seen that tht; maximum volue of are almost same for all 
the values of (~) that are considered, but flattens for 
higher values of ( — ). q increases with increase in (— ) 
and the curve flattens for higher values of (^) [Fig. (3* 10c ) ] . 

A general observation, m free case, car be made is that the 
values of all tho responses are same (or very near that) for 
both the flank widths considered. 

(b) Conical Shells i 

In case of freo conical shells, response w^ increases 
with increase in 9 [Fig. (3.11a)]. At 0 = | wl,57 maximum w is 
equal to 6.65 [not given in figure where as at 0 = 1.31 
it is equal to 0.77 only. Responses Mg and q as shorn 

in Figs. (3.11b), (3. He) and (3. lid) increase with increase 
in 0 and are maximum for © = -g. It can also be seen that 
the curves coincide for the two flank widths considered. 
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'TI 3YI<jHUTRIC PRO J3 Liu i L> 


4 . 1 GENERAL : 

Squat 1 on rj Tor foundation model and shells developed 
m Chapter II have been specialised for thm shallow 
spherical slid Lb and conical snolls oi elastic foundations, 
subjected to antisymmetric loads. 

Shell Lou spherical and conical shells subjected to 
loads m radial direction and moments in radial direction, 
transferred by the column have been solved for fixed, simply 
supported and Croc boundary conditions. The foundation 
thickness has bo^p considered as finite in oil the numerical 
examples that hove been considered. Although, the problems 
with infinite foundation thicknesses, layered elastic founda- 
tions couLcl also be considered without any dxfficulty. 

Results have been presented in non-dimensional form so that 
they can be used lor design purposes . 

In case o C‘ spherical shells two pore-motors which 
accounts fox* thinness of the shell (h/P) and shallowness of 
the shell (R/o) have boon considered. Results have been 
presented for d.i fiorcnl volues of these parameters. In case 
of conical shells results aro presented for different 
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values of Q whicji accounts for the snail ovness of the conical 
shells. 

4.2 EQUATIONS POP THIN SHALLOW SHELnS TO,. ASYl'IHBTKrO 
PROBLEMS % 


4,2.1 S pherical She] Is ; 

In asymmetric problems, the displacement functions 


arc dependent 

on 

both 

r and 

0- 


u s * 

u s 

(r, 

0) 


(4. do) 

v s - 

V„ 

kl 

(r, 

P) 


{4. lb) 

and w _ - 

"o 

( r, 

P) 


(4.1c) 

In vjew of 

the 

Eqns. 

(3.1), 

(3,2), (3,3) and (3.4), 


which arc valid Cor asymmetric case also, the shell equations 
in curv.i 1 unenr coordinates as given in Bqns. (2.25), yields, 


"S u a „ l r u a 

T? + * 7 * 


X. 'll. i £ v a _ x Ha + 2 Hi 
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/ . > r u 3 ift, (1 -i- l s 2) 

+ ( i ■ ii ( J [-j - s -s? ] = — npr“ ~ (4 

r 1 . ?£» . 1 ^ 2u b , x ^ Vq + JL IS 1 + n _ n ) 

*■ "2 £p r &r df> h r 2 ^2 rft £>£ s 

- 1 b V Q 1 -1 1 . 1 C 


r 1 O v o 1 ir , l A. 2-3. + i _ 

x C*2r ^r" ~ ~7 s 2 c>r 2 2r 2? 6 

V e n AW ( 4 “Pg 2 ) 


5 V S lb U S 
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- [ R 
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-- + 


R 


2 c >v s 


4 't 


dir r 


ir + r R + r 2' J J + (1 “^) 
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u 

+ -V + 


isi* i.-: 




. T s i ir r ?_ c^ s 

rR R £ r rR < j3 •* 12 *- 
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2 0 2 ^g ^ 
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0 - - p.--/ 

2 

R h 
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(4.2c) 


where ^/ 2 \ 2 


£4 


- 2 
r 


, 2 c 1 1 £ 4 

< r r~r^ 7vv7 + ?£5 + ?§7 


2 


h 5 


r c r C (3 


■Tffi + 


r 1 , r 

c 4 




r 2 -V^p 2 


4 ' A 

(O' 


r ' c?|S 


(4.5) 


h qua . (4.2a), (4.2b), and (4.2c) represent the equations of 
equilibrium 111 torus oC die placemen! functions u , v and w 

S 3 S 

L'or tli in shallow spherical shell and are applicable to 
m nymme b ric problems . 

4.2.2 Oon icfrl |3 ho).lu : 

I'roccodiu^ on similar lines, as m case of spherical 
sholLo and us uifc Ihjus. (3.1), (3.2), (3.7), (3.8) and (4.1 ), 
the equations of equilibrium, m terms of displacement 
functions, Jor thin shallow conical shell can be written as, 
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r C v 


-L r 

2 - ■ v. 

r o p r 


+ A + Coo 6 ° W£ 


P 


~2~ W S '(1 - H.) [ 


r £ r 
2. 


_i_ Mp , i H 

2 r 2 <* P 2r <jr^p 


_1 

2r 


- <Lla. . Coa & >s , U'O v . , 

2 f> p 2 r >i r ^ S E A ( 4 - 4a ) 

r s 


£ U ( 


- 2 

1 c u 
+ *• 


s2 


G 


f 2 
+ i ..lls . Coa © ffa [1-u ) d v 


' (3 v 2 r> G T p p ^ + 

r • F r C ff ,/ C)P 


s s 

5~T“ 


1 >V o 

^ X j J> 

r ^ v 


*2 

- v + -±_ W . , 1 ^ S I 

2 IJ r 2 £ p ’ r? rc)p ■* 


1 . l^ u 0 -i 


(l-Hg ) 

TT" Y 


(4.4b) 

r U .£ + i- U g . A ^Ts Cos 6 , x ^> u s -, 

L i’ ^ dr 4 r 7TF r w s " (l ~l l s > ~ 3 
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• " - 12 V '. »«, 
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B s F 


r ) o 

who ro V"V is yLvon by J2qn. (4.3) 


Z 

(4.4c) 


4.3 EQUATIONS TOU ELASTIC FOUNDATIONS APl'LI CASES TO 
ASYMMETRIC PROBLEMS: 


Equal ions Cor clastic foundation developed in Chapter II 
aro specialised to suit the foundations for shallow spherical 
sheila and shallow conical shells. 
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In view of the Bqns . (3.11), (5.12), the foundation 
equations in. curvilinear coordinates given hy Bqns. (2.13) 
yield , 
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whe re H is the l"<h ickness of the found etioxi and an me denotes 
derivatives with respect to B is the modulus of elasticity 
of the foundation and p is the Poisson's ratio, Eqns . (4«5 ) 

represents the clastic foundation for asymmetric problems. 

The elastic foundation is characterised by the constants 
which depend ou {y omet no and elastic properties of the body 
[Eqns. (4. 6 ) ] • In this mod ed , eiastic fov lidations of both 
fmi to and inf l m 1v thicknesses and layered elastic founda- 
tions can bo oon..uiorod. 

To!: i ny, onh/ one tor' 1 of the finite series [Bqns * ( 2. 4 ) ] , 
the displacement fu no felons can be written as, 
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In case of externally applied loads on the surface of the body 
(along y = 0), expressions (4.9p) to (4.9r) have to be taken as 
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where P r ( r, p) , P^( r, p) and P ( r, p<) are the externally 
applied surface loads m the radial, circurafeiential and 
normal directions respectively. The body forces P^/ct, p, y), 
Pp(a ? j' y y) find l y( a J P» y) are neglected. 


This gives, 
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Since stresses o , , a , 

*( PY 5 

foundation and flan): and 
it can be seen til*' fc. 



are zero at 
= H cb o os in g 


the surface of the 

j^(y) =X x (y) =<' X (y) 
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4.4 SHELL FOUNDATION INTERACTION E. jUATIOfb i 

4.4 .1 T hin Shallo w S pherical Shell on Elastic Foundation 
S ubjected to A nti-sym ine b rie loa ds ; 

Consider a thin shallow spherical shell on an elastic 
Conndabion carrying asymmetric loads, Ic is assumed that the 
friction between the shell and the foundation is negligible. 
Hence the shell transfers only normal reaction on the 
supporting medium. Therefore, the external load, on the 
shell consists of the known forces Z and the foundation 
reaction f (all referred to unit area). 

z., =, a - p (4.i3) 

As the shell is transferring only normal forces to the 
foundation and assuming the compatibility of the normal 
deflec tions of the shell and the foundation surface, P y 
can be olim mated using Eqn. (4.8c). Considering that 
P (r, (3), PJr, p) are equal to zero in Gqns. (4.8a) and 

(4.8b) and keeping m view the Bqns . (4.12), the resulting 
expressions for the displacements of the shell foundation 
system can be written as [Choosing = 
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where a^, b-^ etc. are given by Eqns. (4.9). 

done id eriny loading to be anti ~f-yw aetnc f the 
displac omont runotions oeu be written as, 
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Substituting Eqno. (-1.15) in Eqna. (4. 1; ) and multiply mg 
suitably with Cor; ( !?n--l ) p cl ,0 or Em. ( 2n •!) (. d ; and integral- 
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where 


a ll J ^11 e ^ Co are given by Eqns . (4.9) 


In all the examples that have beer considered, it is 
assumed that the e3 as tic foundation consists of single elastic 
material and is of finite tlndoiess H. Til- dimensionless 
functions 0^(y), 2 .(k) and ^y) l iav e to b. chosen from the 

physical constrain lx of the problem. to the present inves- 
tigation for finite elastic layer, the '-ollowing function 
are used, 


0 x (y) = ■\ 1 (y) = ‘Rr) = JI i x 


(4.17) 


Substituting 0-^( y ) , Xp(y) end J^(y) fion 4qn. (4.16) in 
Eqn. (4.9) for b^ etc. and integrating from 0 to H , 

the resulting expressions are . 


a 


11 


b 


11 


I 

o 


(1 


""27 

% } 


H 

5 


_ o, 

2(1 + \x Q ) 


(4. 18a) 


1 

H 


(4.18b) 



113 


, _ n. _ 1 

'11 2(1 + |i 0 l 2 


‘11 


j-i n 

0 1 o 


1 


(1 - V) 


11 


H 


2(1 - l - !’ 0 ) 


'11 


f 


11 


S 11 


(1 - u 2 ) 11 
( X'-o> E o H 




27 ? 


ill 


° „ H 

2T 3 


(1 ™ O 


l ll 


& 


I-I 


2U + |lj 5 


k 


11 “ 


2(1 + lO 


1 

2 


ij. E 
“o o 


11 ' (i.. ,r ? y 2 

'O' 


m 


11 


H 

2(1 + Hj 5 


n. 


11 


H 


o a 
H 


2(1 -I- tO 


'11 


2(1 - |ij 5 


( ' .13c ) 

( 4 . lQd ) 

(4. 18e ) 

(4.10f) 
(4.18s) 
( 4 . 18h ) 

(4. 18 1 ) 
(4. 13]) 
(4.18k) 

( 4 . 181 ) 

( 4. 18m) 
(4. 18n) 



115 


To non-- dimen sioiulis* Eqns. ( 4 . 1C ) folLo". m. dime ns ion 
3 e 8 e lu nc t 1 on s e re as a ume d . 
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wher e a is the outer radius oi the shell . .id h is the thick- 
ness ol the shell, 

31 the she 13 is subjected to load u> radia3 direction 
and moments m radio 1 direction transferred by the column., 
then the components ol the force :C, Y and Z are equal to 
aero , 2 . e. 

X = Y = Z = 0 (4.20) 


keeping m view the Eqns. (4.9), (4.19) and (4.20), the 
differential Eqns. (4.16) can be wnttec an fiiven below. 
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Forces and moments m dimensionless form for spherical 
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r dr 


V 


+ (™) (2n-l) -f E1 + (1+H S ) V 
r 


(4.23b) 
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M x (r) 


M 2 (r) 


5 (r) 


M 12 (r) 


q(r) = 


3 h r |i-, o 

--IT - = (fiKle) [- =5 (2n-l) 2 •“ 


n 


+ — 

r dr 


d T_ d 2 7 

n . ii 


+ " T ] 

dr 


(4.23 c) 


Ji'U(r) a , 

-Sr— - (§><!> t 


13 (2n-l) 2 w n + 


dv 


x 


n 


dr 

S(r) 

"Mf" 

dv 


i 2 " 
o T f j- 

+ i» B 72 - J 


(4.23d) 


2(n~7(7j ' b 


(l)(£) 


u 


sn 


v 


+ uxi 

-4? ] 

dr 

r 

M 1? (r)a 
’ TJ 

= (i-i-i.. 

i 

+ jjnzlj 

dw 

] 

dr 

r 

2(1+1^) P. 

rH = 

■ L, 0 


_ ±° i 

_ ^0 


(4.23e) 


h- ■TL ) V 


L 2 n=ll- 

p 2 n 


(4 - 23d?) 


, 1 , ,hwR' r !% c:u : 

Two i? ^ L ~ 2 7^ 


- ( 2n-l ) + 


n 
dr 

a a? 2 


+ ^ 1 (H) !ia 7 i 1 :!? 2 (I) (2„-l)2 

3r a dr L 3r 2 a 


+ 2 <K> *» J 


(4.23g) 
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where 

H 2 (r) 

etc . 

are 

gi veil a s f o 11 0 ’ ■ s 

0 

□ 


N-lCv.P) - i 

yr) 

Coj 

( 2n~l) |3 

(4.24a) 


TT 2 (r,?) = 

I 2 (r) 

Cos 

( 2n~l) p 

(4.24b) 


N x ( r , p ) = i 

I l< r > 

Cos 

(2n -1) p 

(4.24c) 


^ ^2 ( i f P) - i 

■ o ( r ) 

Oos 

( 2ii~l ) , , 

(4. 24d) 


= 

-S(r) 

Ol 1 

( 2n - 1 ) ,3 

(4 . 24 e ) 


H 12 (r,P) - • 

I2 (r) 

Siii 

( 2n-l) 

(4.24f) 

and 

l y (r, p) = 

± y ( r ) 

CoS 

( 2n-l ) J3 

(4.24g) 


(n = 1,2, 

) o • * y 

) 



4.4.2 

Tlun Shallow 

C OjllC 

ol Shells on elastic 

Foundations 


Subjects d to 

Antic 

.vnime 

trie loads ; 



Proceeding on similar lines, as i . > case of spherical 
shell, the yovermnj, differential equations m dimensionless 
form for thin shalLow conic.il shells on elastic foundations 
sub j oo ted bo antis/iuiet rac loads m radial direction and 
^ radial moments tr? ns Terre d by the colu.m ior the foundation 
of single elastic material of finite thickness, can be 
written as given below. 

Differential equations for flank which Is used m 
case of free shell only, are 
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-.2^ 
d u 


fn 


di 1 


n du ^ 1 5 ( 1-Li ) T 0 

J : . —m + [J-, + _11°_ /-■)>- 

r clr 




?r 


( ?:'-i -1) ] u fn •• (-Jn-1) 


(l+p.„) civ. 


2r 


fn 

dr 


( '> '|> ) 

jTisr- <^- x) \'n 


j{Z" ’ 1 ) dw- 

(§) ~~ (4.25a) 

! H dr 


d 2 v 


dr 


fn (2u-l) ( 1+u o ) 0li fa , (2,vl) o } -- 

~ "~~ + "2'"' (ITJ a f n 


v (l-a 0 ) dr 


:r ~ + C “2 + TZ=zry- ~2 

i i d r t x o / p 


„ ?( 2n~l) ^ , 

a \ 2 _ _ JL_ j y 


fn 


>7T^T” (2a-l)(-v) w. 

!(1-|0 h' XXI 


(4.25b) 


d ^ w X-n 

dr 2 


^V 1 ! /a' du ^ 
2 \Y^) ^ dr 


3 V ' : V^' ( a ) “ w fn , 2. ( P.. 1 . (?_) J:a 

2 vr -~\- (- f ) •*-- + 2 ' ~ 

(1 "’V 


i ( 2n-l } (^o" 1 ) ,gn • „ JL ikn. 

+ 2 Ti-v n r ( H } v n - 


r dr 


+ [ + IT-uTj (g) 2 j -'.a (4 - 25o) 


Differential equations for tlie shell -foundrt ion system are 


d 2 u 


dr 


n du 

1 _JL. + [-, 

r dr 


1 + (|)2 + 
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, (1 -V 

‘ 2~r 2 


] u n - ( 2&-1) 


( l+|i Q ) 


2r 


+ ( 2n-l ) 




- 2 

C r 2 r 


v 


n 


\ <>o l) <f> 


cItt 


f 11 


dr 


( - . 26 a) 


,2'~ 
d v 


n 


dr‘ 


tij-ii ^o2 

r O-uJ c1 "' 


1-4- 


(2,1-1)- (3 -V - 

u n 


r v 

JL + 


r dr r 


(3^-D 


r _i 4. j. 2 < an-i) 2 , 

[=2'+ 3(gj + — :=sj 


U-|i 0 )r' 


V 

<s- n 


4 ). -li.Si — - ( 2 n-l) (>) Vi , 


(l-'-' 0 ) 


±:l 


( 4 . 261 ) 


d 2 u 


di 


vSn _ 
-2 


, du -, ( 1-1'v ^ o 

1 „__sn + [1 + — _J 0 S- (2n-l)‘] u 


r o r 
( i+|j. 0 ) 
2 r 


2 r 2 


sn. 


( 2 n-l ) 


dv ( >-|U ) __ 

— _ - . — — — ( 2n**l ) v 
~ + -2 v ' an 

dr 2.c 


a °l- 9 - [H s - ± ] 


V 


sn 


( 4 . 26c ) 


. 2 - 
d v 


sn 


dr 2 



( 2 ji- 1 ) du sn _ 

+ - - — !J -~ 

( 2 n-l) 


r 

ar 


-2 

r 

V 

dv 

sn l 

d r 

[-T? + 

r 

2 ( 2 c-iJ 2 

r 2 U"!\) 

1 "sn 

2 Qps © 

LSSdU. 

w 


T (i 

” V 

-2 

r 

sn 



u 


sn 


( 4 . 26 d) 
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Forces end monento m d omens ion] err fen for conic; 1 
shell con be written cs follows; 


T n ( r ) 

l J ]_( r ) = ““IK 


— -3— ^ (£) [ --sa 

(1 ~ i J s 2 ) ar 


\ X.2n-lJ, u sn , ! n 

+ !L e . \ - 1 — : — ^ '^v, + ~J •' 'O' 


l r 


+ _P - .S-Cos 6-] (4.27a) 

sn r r 


w 2 (r) = 


? r 2 ( r ) 


(k } r ^ 

/- * - sn 

(1 - ) r 


U W 

on , n 


+ God 0 + (j, c 


(■1.27b) 


IL(r)a , a w ( / ~ •, \ 2 

~ / s 1 /li\ r n , I2n-1) 

V r) = -• tb - = (r) [ 7:2" + s V -^ 2 — w. 

dr 1 r 


n r’w 1 
.1 n 1 

-I- - 1 1 


(4.27c) 


r d r 


K 2 (r) - 


:io(r)fl 


(ft) [ . %n->: + i. -J 

r^ 11 r dr 


, 11 i 

+ j. J 

1 dr 


(4.27d) 


3 ,.) _ _ 1 ( h) t .icn-u 3 

[] 12 {r) ~ Eh 2C1HO a L ~ £ 


d v v 1 

sn sn 


(4.27e) 


M 12 (r) 


M 12 a 


, /hv r ( 2n~l) - . ( 2n-l) 

- = (l-n s )(|) [- w n + Hr— 1 - 


(4«27f) 



129 


q(r) = — 


2(H-;. o ) F ;;) 


du 


n 


1 /b x r ' r O a 

Il^i7 V L" ? . 

v 1 r\ ' 


> 1 1 u 
o n 


-M 1 


o 

V 

J! 

r 


dr 


u-o 


0 ( 2 n-l ) ^ (f) 


rA 


2 *> p. 


l a i’ 0 ) ,■- 


1 


(%) { 2n-l ) + 


■ r 
L 


:;( ’'n] 


(4.27g) 


where Jff-^(r), ^(r) ^tc. are .given by JSqnp . (4.24) 
4,5 BOUNDARY CONDI TKVS 


Three types of boundary conditions, namely simply 
suppo rted , Fixed and free types have been considered for 
both spherical and conical shells. 

4 . 5 1 1 B oundary C ondit ion s at the Outer Pounds ry o f Shell ; 

( 1 ) Supported. ; 

The following v r ould be the bounder / conditions at simply 
supported exid 


(a) 

U 1 

= 0 

( 4 . 28 a) 

(h) 

V 1 

= 0 

( 4 . 28 b) 

(c) 

u s 

= 0 

( 4 . 28 c ) 

(d) 

v s 

= 0 

( 4 , 28 d) 

(e) 

V 

= 0 

( 4 . 28 e) 

(f) 

M x 

= 0 

( 4 . 28 f) 
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In dimensionless form, the above mentioned boundary 
conditions for "ootn spherical and conical shells can be 
written as, 


(a) 

u n * 

0 

(b) 

v n 

0 

(c) 

*sn = 

0 

(a) 

v sn 

0 

(e) 

*n = 

0 

/-p \ 

d\ 

L_^J. ~ \X S w n + 


dr r " 




dr 


0 


(4. 29a) 
(4. 29b) 
(4.29c) 
(4.29d) 
(4.29e) 

( 4 . 29 f ) 


( 2 ) Fixed or Built -in : 


The boundary conditions for fixed end are; 


(a) 

^ = 

0 

(4,30a) 

(b) 

\ = 

0 

(4.30b) 

(c) 

u s = 

0 

(4.30c) 

(a) 

v s = 

0 

(4.30d) 

(e) 

W 1 = 

0 

(4.30e) 

(f) 

dw l = 
dr 

0 

(4.3 Of ) 


In dimensionless form* the above mentioned boundary condi- 
tions for both spherical and conical shell can be written as 


(a) 

(b) 

(c) 





0 

0 


(4.31a) 

(4.31b) 

(4.31c) 



0 
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(d) 

v =0 

sn 

(4.31d) 

(e) 

w =0 

n 

dw 

(4.31e) 

(f ) 

™ = 0 
dr 

(4 . 31 f ) 

( 3 ) Free Case s 


Boundary conditions at the far end of the 

flank can 

he written 

as 


(a) 

% = 0 

(4.32a) 

(b) 

v^ = 0 

(4.32b) 

(c) 

n = 0 

(4.32c) 

Boundary conditions at the free end of the shell 

written as. 

can be 

(a) 

Uq_ - u f 

(4 33a) 

(b) 

V 1 * v f 

(4.33b) 

(c ) 

W 1 ” w f 

(4.33c) 

(d) 

Normal force m flank 



= Notmal force m foundation 

(4.33d) 

(e) 

Shear force inf flank 



= Shear force m foundation 

(4.33e) 

(f) 

Noimal force m shell = 0 

( 4 . 33 f ) 

(g) 

Shear force in shell = 0 

(4.33g) 

Ox) 

Transverse shear m shell 



4- Transverse shear in foundation 
= Transverse shear in flank 


(4.33h) 



r 


(l) M x = 0 


(4.33-0 


Above mentioned boundary conditions m dimensionless 
form for sjiierical shell can be written as given below. 
Boundary conditions at the far end of the flank are 


(a) 

(b) 


u 


fn 


fn 


( c ) w, 


fn 


= 0 
= 0 

= 0 


Boundary conditions at the free end are 


(4.34a) 

(4.34b) 

(4.34c) 


(s ) 

“n = 

ll fn 

(4.35 a) 

(b) 

^n = 

v fn 

(4.35b) 

(c) 

w n = 

w fn 

(4.35c) 

(cl) 

du„ 

fn 

dr 

du u _p u / „ , * 

__,n ... r _f n + L2ji=ll ^ 

» + P 0 L _ _ + -- ^fn 

dr r r r 



- ^ + | (|) vr fn - | {§) w n ] = 0 

(4.35 d ) 


dv P dv 
fn n 


(e) - L&rp- ^ fn + u n + - 

r 1 r dr dr 


V, V 

Jk + = o 


(4.35e) 


(f) 


du 


'sn 


dr 


(|) h ♦ ^ \, + a -f + <l> = 0 


(4.35f ) 
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dv v , 

<e) -12. --S5L _ iziizll - 

dr ? v sn 


= 0 


(4.35s) 


(h) 


(i) 


6(1-1^") L dw, 

e“ 3 c 2 U fn - I % + J <f> 


(1 + l‘J 


dr 


2— 

1 W v, -- o ~ d"V , d W 

1 - n n /lio/io r ni , 1 n 


f <§> + (S) , (f) J [-3^+ i 

d x c.' r r 


3 


n dl7 

A a 

•■2 i~ 

r dr 


dr" 


( ^n-1) 2 du n , .2 / o n n n 2 “ -i n 

rr + ( 2n-i) w j = o 

d r r J n 

( 4 . 35 h ) 


d \i u dw 

■n I cj 

+ 


n ‘ ' S rr- - (2n-l) 2 w„ = 0 (4.35i) 


dr 


dr 


n 


Boundary conditions m free case for conical shell m 
dimensionless form can be written as given below. 

At the far end of the flank the b oundary c auditions 

are „ 


(a) u Xn 

=r 0 

(4.36a) 

< b) v fn 

= 0 

(4.36b) 

(°) w rn 

= ao 

( 4 -36c ) 

Boundary 

condi bions at the 

free end of the shell are „ 

(a) ^ 

rf 

tH 

II 

(4.37a) 

(M v n 

= " V fn 

(4.37b) 

(c) w n 

= w fn 

(4,37c) 
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du, 

(d) -ia 


dr dr 


dr r 

( 2n-l) 

v fn 

\x 

S J 

~r 

i l o - 

!- l 0 _ 

3|’ 0 

. O V 

r u fn “ 

~r u n + 

”2 

<ij) 


(2n-l) v. 


n 


2 l^o ( H } = 0 

U, 


(A . 37d ) 


(e) -( 2n~l ) + ( 2n-l ) ^ + fZjfc _ fn _ l£n 

r r dr dr r 


v 

+ * 7 T = 0 

r 


(4.37e) 


du f 

dr 


v sn u sn „ T 'n 


(£) -^ + |j [(2n-l) -fii + _Sii + _a cos 9] = 0 (4.37f) 


(g) 

-( 2n-l ) + 

r 

dv 

sn 

dr 

(h) 

12dV> E o 

2(1 + \i Q ) B s 

(§;)3 

l h' 


- i (^) 

3 V 

d ”fn 

dr 


n \2 dw 

■ A jljcrjr,.) — \ 

— 2 , — 

r dr 

(i) 

d 2 w u dw V 

n . s n ' 

dr 2 r dr 



= 0 


r 

u u. 


(4.37g) 


] - 


dw 

dr 


a \ 1 e \ , ifTa 

dr^ r dr 2 r 2 dr 


w n = 0 (4.37h) 


w u = 0 


(4.371) 



4.5.2 Boundary Go ndit ions at the Column E nd s 

The antisymmetric load distribution along the ctiLumn 
has been taken 

°l> = ^ 7 tF' Oorj ( 2n-l) j3 (4.38) 

where P is either radial load or moment in radial direction 
and b i-> the radian of the column. If the column is assumed 
to be very stiff or rigid m com pans an to shell, the load 
distribution along the column can be token as, 

% = ib 308 P (4.39) 

since no higher order mode shapes are possible. 

Jj oundary conditions at the column end can be written 

m three ways. Firstly by assuming that the foundation core 

(por non of foundation below the column) is rigid or very 

stiff. Secondly by assuming that the core of foundation is 

{ 

of negligible stiffness and lastly by assuming a linear 
combination of the above two extremes. To be on safer side 
second criterion has been considered m writing the boundary 
conditions . 

4.5.2. 1 Load m Radial Directions 

Boundary conditions at the column, end of the shell 
for loads m radial direction can be written as? 



(a) formal Cores in shell ( If ^ ) 

=■ Load per unit length 
dw 

(b) - o 

dr 

(c) Transverse shear in shell 

+ Trai averse shear m foundation = 

(d) Normal force m foundation = 0 


(4 .40a) 
(4.40b) 


(4. 40c) 
(A .40d) 


(e) Shear force in Shell = 0 


(4.40e ) 


(f) Shear force in foundation = 0 


(4.40f ) 


In dimensionless fom the above no ntioned boundary 
conditions for spherical shell can be written as, 


(a) + (f) * a + „ v sn + hn + (|) v )] 

dr r r 

Cos (2n-l) J0 = -»- ■— (#) — ^ G os p (4.41a) 

TC Op, h 

s 


from orthogonal properties of cosine function, it 


can be shown that, 


sn , /a 


+ ( jf ) w n + p s [ ( 2n-l ) 


v u _ 

sn , sn — 1 

— + + ^ V 


71 D E li 


n = 1 


0 


n > 2 


(4.41b) 
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( 1 ) 


dv 


n 


0 


c*r 


(c) <£) 3 (f-K [ 


'R 


d^v 


n + 

dr' 5 r dr 2 


dw 


u 


-2 i“ 

i’ d r 


(4.41c) 


aw 


( 2n~l)“ ~ v, n L 3 , r . ,2 - , 
' i — h -?■ (^xi~l) v^] 


-2 


6r 


12(1-^ ) B 
E 


. s -- o r In , 1 /ll 

2(1 + uj E s 2 


-| TT dVJ 

+ 9 (?) T 2 ] 


3 'a' 


dr 


= 0 (4.41d) 


(d) 


du„ 


n 


- + 


dr 


!' ( 2n -l) u 4n 

? n + 1*0 ^ + -r° <S> *» 

r r 


= 0 (4.41e) 


(e) 


n \ _ dv 

u + — - 

' ■* * % ^ 

dr 


v 


yn 


_sn 

r 


- = 0 


(4.41f ) 


(f) 


( 2n-l ) ~ , 

i L u + 


dv 


n n 


n 


= 0 


dr 


(4.4 lg) 


For conical sheila, boundary conditions in dimensionless 
form can be written as. 


, , du sn „ u sn 

(a) —z- + |i -zr 

dr r 


v. 


+ (2n-l) |i c 421 + , is cos © 


w 

r 


^“O 4w P 


= 0 


n = 1 

% 0 E 

s 

n h> 2 


(4.42a) 


(To) 


dw. 

dr 


n 


= 0 


(4.42b) 



(e) 


3- 2-- - 

d w d t > j dw /r», n \ 2 dw 

— Ji _ JL _ — ja + Jk a 

dr r dr^ r dr r dr 


. 2L2n=lii - + !o ( a )3r % 

-3 'n + IThi 0 ) L s V I- 2 


H Tj dw 

+ 1 rt) Hi j =0 

' tJ dr 


(4.42c ) 


du v„ n, - „ _ . 

(d) ~ + r. 0 (2n-l) ---+ -S u + |l'- (w) w =0 (4.4 2d) 

dr 0 r - n 2 o a n 


(e) - ^ 


_C2n- ; lJ ~ + 

r ' sn 


dv. 


sn sn 


dr 


= 0 


( f ) „ iinJL) 


dv v 


r 


u n + 


n 


d r 


n 

r 


= 0 


( 4. 4 2e ) 

( 4 . 4 2 f ) 


4,5.2. 2 Moment m Radial Directions 

Here again only first mode is effective „ .Except for 
boundary conditions (a) and (b) of Bqns „ (4.4 0) all remaining 
four boundary conditions are same as m case of load m radial 
direction . 

These two boundary conditions ore, 

(a) Mi = M , . , m per unit length (4.43a) 

v ' 1 applied £ 

(b) Normal foice in shell = 0 (4.43b) 

For spherical shell, the two boundary conditions in 
dimensionless form can bo written as 
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, 2 ~ 

p u cl VJ 

(•) c-=^ + 1*. 

dr k 


■%-U 2 ;; + i faD] 

—d A — ~ , J 

r r dr I 


'' 7t ^ b^ D n ” 1 


= 0 


du r 

dr 


n > 0 

iSS=il v 


(4.44a) 


(b) _3 + (|) ; a + ^ -V S11 + + (g) T , n \ ] 


= 0 


(4.44b) 


where M is the applied moment 9 

D is the rigid it j. modulus. 

(^■)(^) has been taken as the load parameters. 

Similarly for conical shells these two boundary 
conditions are 

dr | r r dr 

«-£[<!><§>] i"- 1 


n 


= 0 


jn = 2 


(4.45a) 


W 


~ + i*„ — 1 


V _ XT 

i A . n A 


dr 


t (ai-1) u CJos 9 = 0 

r e r r 

(4.45t>) 


4.6 NUMERICAL CALCULATIONS AND CONG HJSIOFS : 

In all the numerical examples that have been solved, 
following parameter values have been kept constant. 



§ = 2.00 


b 

a 


= 0.2 


= 0.25 


|i 0 = 0.17 


R, 


In case of spherical shell, (— ) = 4.00 has been kept 
h Ec 

constant while (£) ?nd (—■) have been varied. In second set, 

R L ° v 3 S 

(") = 0.0333 has been kept constant while ( j) and (g— ) have 

o 

boen varied. Similarly m case of conical shells (^) = 14.0 

B 

is a constant while 9 and (gr-) take different values. Results 


have been presented in plots drawn between r and functional 
values at unit load parameters. load parameter m case of 
radial load is taken as, 


= <f>< 


S h 
s 


■) 


(4.46 ) 


and in case of moment is taken as 


i = (fug) 

Pro in boundary conditions at column end it Can be 
concluded that only n = 1 gives rise to response m struc- 
ture. For n >2, there would not be any loading on structure, 
hence there would be no response of the structure. 
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.Keeping this in nine!, response functions per unit 
value of lo ad parameter have been defined as [writing n = 1] 


U 1 

= 

•g/r 

(4.48a) 

U S 


u^/F 

(4.48b) 

w^ 

•= 

l Out 

< 

1 1 

IIS 

(4.48c) 

% 

= 

fg/p 

(4 .43d ) 


- 

? 2 /p 

(4.48e) 

M i 


1 — 1 

(4.48f) 

«2 

- 

ii 2 /p 

(4.4Qg) 

q. 

= 

i/p 

(4.48h) 

Load 

in 

Radial Direc tion i 



4 .6.1.1 Spherical Shell; 

(a) Simply Supported and fixed Boundary Conditions; 

As shown m Figs, (4.1a), (4.1b), (4.2a) and (4.2b) 
responses u g , increase with decrease m (^; also maximum 
value of responses u g , w^ in fixed case arc slightly lower 
than those in simply supported case. iT^ remains positive all 
along r and increases with increase in (g) as shown in 
Figs. (4.1c) and (4.2c). M for simply supported case is 
shown in Figs. (4. Id) whereas is shown m Figs. (4.1e). 
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Af shown in Fi;s. ( 4 . 3 a), (4.4a) vt lues of responses u g 

"D 

increase with decrease m (~) . w, decree ses with increase 

El X 

FJ 

m (-j|) and the curve Clottens for higher values of (~) 
which shows the tendency of v going to zero along r for 
large value of (^), that -is, in case of 0 plate [ Fig. ( 4 . 3 b ) , 
(4.41,)]. Nj remains posit I've for all values of r and 

p 

increases with increase in (~) and tends bo become maximum 

a 

for a plate [Fig. ( 4 . 3 c), ( 4 . 4 c)]. From ligs. ( 4 . 3 d) and 
(4.4c), lb can be seen that 'q decreases wi th increase of 
(~) and it tends to zero along r as (~) takes higher values. 
Similar conclusions can be dra wn from Fi,l . ( 4 - 4 d ) showing 

variation of 11^ along r. 

(b) Free Boundary Condition; 

is zero at the far end of the 1 lanL and it is 
negative near the free boundary and becomes positive away 
from the free bound* iy which would obviously oe the result 
in case of free she Us [Figs. ( 4 . 5 b), ( 4 . 6 a)]. w approaches 
zero along r for larger values of (~) ^ can be seen xrem 
Fig. (4.6a). *q is slightly negative near the free bounds, ry 
and positive in the remaining portion [Fi/e. ( 4 . 6 b)]. 

It can also be observed that q tends to zero along r for 

large value of (~). 
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4.6,1, 2 Conical Shell: 

Simply Supported and Fixed Boundary Condii ions : 

In Figs, (4.7o)» (4.8a) variations of u. & have heen 
shown for different: values of 0 for sum ply supported and 
fixed boundary conditions. Variation of 1'-^ has been shown 

■'"'w 

m Figs. (4.7b) and (4.8b) for different values of 9. 
remains positive along r and is maximum at the column end 
where it is equal to the prescribed IT-^, H-p as shown in 
Figs, (4.7c) and (4.8c), is zero along r when 0 = (plate), 
which would obviously be the case since no moment can 
develope if plate is subjected to inplane forces. 

4,6.2 Moment m Radial Direc tion ; 

4.6,2. 1 Spherical Shell 

(a) Simply Supported and Fixed Boundary Conditions: 

< w , -^ ? as shown m Fig. (4.10a), rerns inc positive all 
along r and the value of w^ increases wi bh decrease of 
(~) ratio for fixed case. Similar is the Ccse for simply 
supported boundary condition (Figure is not shorn). Varia- 
tion of MU along r is almost same for both simply supported 
and fixed boundary conditions and for different (jj) values, 
except near the boundary where M r is zero in case of simply 
supported shell and is having slight positive value in 

is maximum near the column end where its 


fixed core. 


value is e qua 1 to the 


14 *. 

prescribed moment [Figs . (A . 9 ) , (4 . 10b) ] . 
wj_ is shorn m Pi,'. (4.12a) for riff ©rent vs lies of 

p . < ^ — 

(~). li-j is ogam nrving almost the sou* vc riation along r 

■p 

tor all values of (-„ L ) consider’d and for both simply supported 

Ct 

and fixed types of boundary conditions (except near the 
boundary)* Only li^ for fixed boundary condition has been 
shown m Fig, (4.1-^b). Mg for simply suwoorted case has been 

*^J 

shown m Pig. (4.11). It can be seen that fig increases with 
increase in (--■) ratio. 

C*i 

(b ) Free Boundary Condition: 

\ J ^ 

Fig. (4.13a) shows the variation of alone, r for 
different values of (p) . jt can be seen that increases 
slowly from zero at the far end of the f lei ft to borne small 
positive value at free boundary of the slie 11 and goes 
negative away from the boundary. Fig. (4 .id a) snows the 
variation of wp along r for different vain es of (t^). Similar 
observations, as Jr .i been made above, can also be made. 
Variation of along r is same (or very near that) for all 
(».) [Fig. (4.13b)] variation of Mg along' r for different (— ) 

J L 

is shown in Fig. (4.14b). It can be observed that Mg increases 
with increase in ( — ). 


4.6.2. 2 Conical iJliulls 


14 5 


(a) Simply supported and fixed boundary ccr it Ions ; 

Variation of w^ along r lias been s’io T, r in Figs. (4.15a) 

* , 'S^ 

Find (4-16) for dilfeient value.? of 6. 1S maximum for 
© = ■- (plate). As shown in Figs. (4.1^’‘ 1 ), li^ is almost 
same for all values of 0 considered. Similar is the case 
for fixed boundary condition (Figure is not shorn), M-j_ 
attains maximum value at the column end where it is equal 

rN/ 

to the prescribed value of 

(b) Free boundary condition; 

Fig. (4.17) shows the variation of w^_ r for 

different values of 0. It can be observed from the plot that 
is zero at the Car end of the flank and slowly increases 
to some positive value near the free boundary and remains 
positive over a lai go portion of shell and rap idly drops 
to a negative value near the column end. lb is negative 

71 

value is much larger than its positive value. At Q = 
w has maximum negative value. 
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CONCLUSIONS AND ESCOilMBJ DAT IONS FOR FURTHER WORK 


5.1 COfTCLUSIONSs 

Foundation nod el in orthogonal cu’ulinear coordinate 
system, m which the displacements m all the three ortho 
gonal dire o t ions have been included, has eon developed 
using vana to o nol method, A wide vsnctp of problems of 
clastic foundatiois can be solved using this model because 
oL‘ its general nature, 

j roblcms oL shallow spherical and conical shell on 
clastic foundation;, have been solved. AJ 1 the three 
important typos oJL boundary conditions, namely simply 
suppor tod , X’j^cod and free, hove* been co n f id one d. In case 
of free boundary conditions flank (found: Lion beyond shell) 
ho vo bo (Mi conoid om’d and on effecti-ve Clan 1 ' width which 
needs to bo considered has been suggested. 

JvoJiuns mu Lin segment technique when has been used 
to solve the boundary value problems, has worked very well 
sud results obtain'd during forward integration have matched 
upto so von rj i gmf J can t fjgn.ro the results obtained on bac to aid 
ml, igrn tion. 
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The result:. 0 obi.:. ined for the numeric r 1 exam lee have 
been discussed at th^ end of the Chapter III and TV. Fumeri- 
cal results for all the variables and for all the cases 
considered have not been presented 7 howevn important results 
have been presented to jive both, quant its 1 ive and qualitative 
presonta bion of forces and displacements tor oil the cases. 

Although no cost analysis has been done but it can 
bo easily seen from J'Lps. ( >.7), (3.8) and (3.0) that shells 
would be much more efficient as foundation of column than 
the plate,; are . 


One important observation that c-'n he made from 
(3.3a). (3.3b), (3.4a) and (3.4b) is that the shell 

Q lion over the foundation. This happens uo cause no friction 
has boon assumed between shell and foundation. It can be 
noted from the variation of q that it is not constant over 
the depth of Lcnmdnt ion. 


In the numriicnl - xrmplos , tliicto'oss of the founda- 
tion has been considered as finite and the foundation hcis 
been assumed to be made of nxa^lc material but foundation 
ol InlltiiU: ihsoWnri and of different clastic layers can 
be,- oounulorud with equal ease by proper choice of functions 


of normal distribution of displacements. 
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5.2 KSCOiiIUo'DAI 10 L' v'URIjIhR i/ORi; r 

In the c > all o 1 s vuvv follovm^ s to 'eh 3 s re ss where 
ieiu.aroa vcmV ere -'nuircd. 

( z ) Thiele sir 1J; or classic foundation s . 

( , ) Dynamic ? .missis ol shells on classic foundation, 

(c) Us i ol' Hid proposed model on othc i tyocs of sheila 
J i • Uypar, 

(d) fo consider cn’tiun cuounb of Criction between slid 1 
i uio found r bum to make the anaLysi snore realistic. 

In nut shoJl, ,'io or] amount o l rerjooro 1 work , both 
f.violytjc 1 end exp i'oonb'1 is required on shells on el so tic 
found t bun'.f btfor iloi,i(jn charts can be prop? red. 
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